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The Behaviour of the Cathode Spot on an Undisturbed 
Liquid Surface of Low Work Function 


By K. D. FROOME* 
Department of Physics, Imperial College, London 


MS. received 19th September 1949 


ABSTRACT. The cathode spot of transient arcs using the liquid sodium—potassium alloy 
for cathode has been studied by means of the Kerr cell apparatus capable of taking a 
sequence of photographs. It is found that for a given current the spot takes, as in the case 
of mercury, the form of a line or broken line of total length roughly proportional to the 
current. This is also true for a spot in a magnetic field. Again as for mercury, if the 
current rises slowly the length of line increases proportionally to the current. If the rate of 
growth of current is greater than a value lying approximately between 10° and 2 x 108 
amp/sec., then fresh spots form and spread out radially from the point of formation into thin 
semicircular or circular lines moving with a maximum radial velocity of about 104 cm/sec. 
The apparent current density of emission lies between 2°5 and 5x 10% amp/cm?, and 
consequently is higher than for mercury. 


§1. INTRODUCTION 

HE results described in this paper have been taken by means of the Kerr 

cell apparatus capable of taking a series of superimposed photographs 

which show a number of positions of the cathode spot as it moves over the 
cathode. The interval between each exposure can be varied from a fraction of 
a microsecond to many microseconds, and the number of exposures of the 
cathode spot taken during a discharge can be adjusted from one to as many 
as are required to build up a complete picture of the behaviour of the emitting 
area. The exposure can be varied independently of either interval or number. 
This apparatus has been described elsewhere (Froome 1948). 

Observations are made by means of a low-power microscope arranged to 
look down on to the liquid cathode surface, the Kerr cell itself being placed 
just outside the microscope eyepiece. ‘The actual optical arrangement used 
is shown in Figure 1 of a paper describing details for an undisturbed mercury 
cathode (Froome 1949); reference should be made to this paper for other 
allusions to mercury cathodes. 

The transient discharges studied are obtained by discharging a single 
condenser, or a number of condensers, in the form of artificial lines, charged 
to approximately 100 volts through the tube and series impedance. ‘The former 
is used when varying currents are required, and the artificial lines can be used 
to produce constant current pulses of 10 to 450 amperes. 

* Now at the National Physical Laboratory, Teddington, Middlesex, 
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The discharge is termed an arc because the cathode spot emits a large number 
of electrons, and because there is a low overall potential drop. The term cathode 
spot is used to denote the whole emitting area on the cathode irrespective of its. 
exact form. 

The liquid alloy used has a constitution of approximately 25%, sodium and 
75°%/, potassium, but almost identical results can be obtained with any other 
proportion provided it is liquid at room temperature. The alloy has been chosen 
for two reasons. Firstly, it is liquid at room temperature, and a liquid surface 
is desirable because it is reproducible; secondly, it has a low work function, 
whereas mercury has a high work function (4-5 volts). The curves given by 
Thomson and Thomson (1928) for photoelectric emission from the alloy indicate 
a work function of about 2 volts. 

In order to use this alloy, the discharge tube must be sealed off. It was not 
found possible to purify the alloy by the technique of evaporation usual for the 
alkali metals, for this separated the constituents of the alloy and deposited 
potassium in an irremovable film over the discharge tube window. The following 
technique eventually proved entirely satisfactory. 

The tube is made with a long thin glass pumping line attached. ‘This line 
consists of a number of bulbs of about one inch diameter separated from each 
other by a narrow capillary neck. The end farthest from the discharge tube 
(and nearest the liquid air trap) is fitted with a cone joint device by means of which 
the previously prepared and consequently rather dirty alloy can be introduced 
into the vacuum line. The whole is out-gassed as thoroughly as possible and 
filled with argon to the operating pressure (0-1 to 0-01 mm.). Then the liquid 
alloy is tipped into the bulb farthest from the discharge tube and the whole 
sealed off from the pumps. ‘The alloy is then passed from one bulb to the next 
through the narrow joining capillary tubes, and by the time it reaches the discharge 
tube it appears brilliantly clean. The discharge tube is then sealed off from the 
cumbersome remaining line of bulbs. 

The type of discharge tube used is almost the same as that described earlier 
(Froome 1949, Figure 1), but no amalgamated copper foil is wrapped round the 
inside of the cathode neck. Consequently the alloy surface is always somewhat 
curved, and for this reason some of the photographs shown are in focus only in 
one or two successive exposure positions. 

Finally, the cathode is cleaned further by running the tube for several hundred 
flashes at high peak current. ‘This is absolutely essential if reproducible results. 
are to be obtained. The cathode spot has a strong cleansing action upon the 
cathode. 

In all the Kerr cell photographs reproduced the cathode spot has started 
from the bottom of the picture. Hence each photograph shows a number of 
positions of the emitting cathode area separated by the time interval stated, with 
time increasing as the cathode spot moves upwards and, in the case of rapidly 
rising currents, outwards as well. ‘The discharge is initiated by means of a brief 
high tension pulse applied outside the discharge tube near the cathode. It is. 
found that the cathode spot almost invariably forms on the liquid surface adjacent 
the glass wall of the tube nearest the position of application of the triggering 
pulse. In some of the photographs (notably those shown in Plate II) the curved 
surface of the glass cathode wall can be seen in black outline along the bottom 
of the picture. 
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§2. RESULTS WITH CONSTANT CURRENT PULSES 


It is difficult with this alloy to obtain results from an undisturbed surface. 
‘The emitting area rapidly breaks up into a large number of minute separate 
spots. ‘These ‘unit spots’ seem to form in self-made dents in the surface. The 
same behaviour is observed for mercury, but takes considerably longer. 

Figure 1 (Plate I) shows a photograph with exposure 1 microsecond taken 
4 microseconds from the start of a constant current discharge of 410 amperes. 
It is seen that even in this short time unit spots have formed, whereas the 
corresponding time for mercury would have been 20-30 microseconds or more. 

Such observations suffer from the limitations imposed by these dents (Froome 
1949); and much of the work for constant current pulses has therefore been 
undertaken with the cathode spot in a strong magnetic field, since this produces 
a movement of the cathode spot sufficiently rapid to inhibit unit spot formation. 

For the cathode spot in a magnetic field we again find, as for mercury, a 
marked reluctance to travel with a velocity of more than 104 cm/sec. Also, if 
the velocity is much below this value, excessive formation of unit spots takes 
place, and these are left behind the main emitting area until they eventually 
extinguish. 

Nevertheless, there is a marked similarity in behaviour to that for mercury. 
Figure 2 (a)—(d) (Plate I) shows sequences of exposures of the cathode spot 
for a constant current of 90 amperes at two magnetic field strengths and with 
two different exposures. ‘The magnetic field is across the picture, the spot 
moving in exactly the opposite direction to the ponderomotive force in the 
emitted electrons. In Figure 2(a) the field is 900 gauss, in (4) it is 1,800 gauss. 
The exposure in each case is 0-1 microsecond, and the interval between exposures 
is 4-2 microseconds. At the higher field it is seen that the emitting area is in 
the form of a line similar to that for mercury, but even at this field strength, 
with the consequent high velocity of the spot, there is evidence of the formation 
of a few unit spots which have been left behind in the wake of the main emitting 
area. At the lower field the line has become extremely ragged, or broken, due 
to such spot formation, and thus it is dificult to assess the actual emitting area. 

Figure 2(c) shows a sequence, with longer exposures (0:3 microsecond) but 
with the same interval between each, of another discharge again in a field of 
900 gauss. The raggedness is more clearly seen—as are the numbers of unit 
spots. Figure 2(d) shows a spot in the higher field (1,800 gauss) with this longer 
exposure. The line-like structure of the main emitting area is again seen, but 
there is also evidence of unit spots which have broken off and been left behind. 
This behaviour is due to the dents such unit spots apparently make in the liquid 
surface. Once a dent has formed the spot is not capable of rapid movement ; this. 
applies to mercury also. 

Figure 3 (a)-(c) (Plate I) shows multiple exposures of cathode spots carrying 
various currents in a magnetic field of 1,300 gauss. ‘Ihe exposure in each is 
0-3 microsecond, and the interval between consecutive exposures is 4:2 micro- 
seconds; (a) shows the spot for a current of 14 amperes, (5) for 60 amperes, 
and (c) for 165 amperes. The emitting area is seen to be chiefly a line for the 
two lower currents, but for the highest it consists of a very ragged (broken) line 
plus a large number of unit spots in its wake. It is also seen that these spots 
have a life-time of the order of the interval between consecutive exposures,. 


namely about 4 microseconds. 
2 D-2 
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Figure 4 shows graphically a plot of the current against length of this emitting 
line. As for mercury, this curve is obtained from a large number of different 
results. In this case, however, only those results which show the best line- 
structure have been included, and a correction has been attempted for the 


probable effect of the residual unit spots. 
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As for mercury, we see that the length, /, of the emitting line is proportional 
to the current, 7 (amp.), where /=Kzi cm., and K=4-1 10-4 cm/amp™. This 
is under one-half the value obtained for mercury. 


§3. RESULTS WITH VARYING CURRENTS 
The results with varying currents show a strong similarity to the corresponding 
behaviour of mercury. As before, these results are obtained with the spot not 
in a magnetic field, in order to maintain simple conditions. 
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Figure 5 (Plate IT) shows the appearance of the cathode spot for the electrical 
conditions shown in Figure 6. It will be remembered that for mercury these 
conditions produced the formation of fresh spots. Figure 5 shows three exposures 
of 0:3 microsecond separated by 4:2 microseconds, the first being taken 2 micro- 
seconds from the initiation of the discharge. It is seen that only the first shows 
anything of a line-like structure, and that in the second and third slow-moving 
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unit spots have formed; consequently in each of the last two exposures one can 
scarcely distinguish the positions of individual spots. 

We have already seen that for a spot in a magnetic field the length of the 
emitting line for a given current is under one-half that of mercury. Hence, 
since the spot is reluctant to move much faster than 104 cm/sec., one can expect 
that at very high values of rate of growth of current the cathode spot will be 
forced to move in the form of a semicircular line at its maximum velocity, or 
to form new spots. 

When the emitting area does take this form of a semicircular line moving 
outward with a maximum permitted velocity of 104 cm/sec., its maximum 
length L after t seconds will be given by L=1047t cm., or, after ¢t sec., the 
maximum current 7 emitted by such a line will be given by L/K where K 4 x 10-4. 
Hence the maximum current 7 =8 x 107 t, and therefore di/dt =8 x 10’ amp/sec. 
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represents the maximum rate of growth of current which we can expect to be 
maintained by a single semicircular spot. (The spot is almost invariably 
semicircular since it starts adjacent to the wall of the glass tube and can therefore 
spread out only into the tube. If the spot does form on the cathode surface away 
from the tube wall, then it can form a completely circular line, so that the 
maximum rate of growth of current sustainable by such a spot will be twice the 
above value.) 

Figure 7 (Plate II) shows a sequence of four exposures (exposure 0-1 micro- 
second, interval 1-0 microsecond, the first exposure occurring just under a 
microsecond from the start of the arc) of the cathode spot obtained by discharging 
a 20 wr. condenser initially at 60 volts through the tube in series with a total 
inductance of about 0-3 ~H. with negligible total circuit resistance. ‘The peak 
value of rate of growth of current is 108 amp/sec. The cathode spot has formed 
as a single spot adjacent to the tube wall at the bottom of the picture, and then 
disintegrated into a semicircular line moving radially outward with a velocity 
of about 104 cm/sec. The characteristics and time of the discharge are shown 


in Figure 8. 
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Figure 9 (Plate II) shows the result with the same exposure conditions, with 
the condenser charged to 108 volts. The peak value of rate of growth of current 
is now 2:4 108 amp/sec. The actual values of current and rate of growth 
of current can be obtained from Figure 8 by multiplying by a factor 2-4. The 
time of the discharge is still as shown in Figure 8. In Figure 9 two spots have 
formed, disintegrating into semicircular lines moving with a velocity of nearly 
104 cm/sec. Five positions are shown, the first appearing about 0-5 microsecond 
after the start of the discharge, so that just over 4 microseconds of discharge 
time are shown. 

Figure 10(a) and (b) (Plate II) shows two photographs of spots obtained 
with the same exposure and circuit conditions as before, but with the condenser 
charged to 144 volts. The time of the discharge is still as shown in Figure 8, 
and actual values of current and rate of growth of current can be obtained by 
multiplying the values here shown by a factor of about 3-4. In (a) two spots 
have simultaneously formed and spread out into lines moving with an_ initial 
velocity of about 1-3 x 104 cm/sec. The cathode spot is shown for about the 
first five microseconds in about six positions—the first exposure of the right-hand 
spot occurring practically simultaneously with the formation of the spot, the 
left-hand spot having formed about 0-1—0-2 microsecond later. 

In (b), three spots have formed simultaneously and spread, initially moving 
outwards at 104 cm/sec. ‘The first few microseconds are shown in three or four 
exposures, the first exposure being about 0-1—0-2 microsecond after the start. 

Hence we see that a spot can support a maximum rate of growth of current 
of about 10® amp/sec. (if it starts outwards from the wall of the tube). If higher 
rates of growth of current are required, then extra spots appear simultaneously. 
This is approximately what was predicted from the study of the spot in the 
magnetic field. If the spot should form away from the wall of the tube then the 
critical value would be 2 x 108 amp/sec. 

The formation of fresh spots should not be confused with the appearance 
of large numbers of unit spots formed when the velocity of the emitting line 
decreases to slightly less than 104 cm/sec. ‘These unit spots are not fresh spots, 
but broken-up pieces of emitting line. Fresh spots only form under the conditions 
of rapid growth of current. If the current is relatively slowly rising (as was 
the case in Figure 5) and unit spots have formed, the size of the emitting area 
increases by the breaking up of unit spots into greater numbers, and not by the 
spontaneous formation of fresh spots. If the current is relatively slowly rising, 
but the velocity of the emitting area is sufficiently high for unit spot formation 
not to take place (e.g. with the spot in a magnetic field), then the emitting area 
preserves its line-like structure, increasing proportionally with increasing current. 
This suggests that when unit spot formation has taken place each-unit spot 
consists of a piece of line in the bottom of a dent, and that if the current is increasing 
the pieces of line in some of the dents increase, resulting in instability of the 
dents, so that they break into two or more separate dents. If the current is 
decreasing the total length of emitting line decreases proportionally, and the 
number of unit spots decreases. 

Returning to the discussion of rapidly rising currents with the cathode spot 
not in a magnetic field, Figure 11 shows a plot of length of emitting line against 
current for a number of discharges similar to those shown in Figures 7 to 10 
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(PlateIT). Again we find approximate linearity between the length of line, /, and 
current, 7 (amp.), such that /=Ki cm., where K =3-7 x 10-4 cm/amp"t. 

In view of the experimental difficulties, this value of K may be taken as 
identical with that obtained for the spot in a magnetic field. 


§4. DISCUSSION AND CONCLUSIONS 

It is a curious feature that, whereas for mercury (Froome 1949) the fresh 
spots form in front of the previous ones at high values of rate of growth of current, 
with the liquid alloy they almost invariably form along the edge of the cathode 
adjacent to the glass wall. Further, there is less time lag between the achievement 
of the high value of rate of growth of current and the formation of the fresh spots 
required to sustain it. For mercury this time lag may be of the order of a micro- 
second, whereas for the liquid alloy it appears to be nearer 0-1 microsecond. 


[ 
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Figure 11. 


Most of the conclusions given for mercury apply also to the liquid alloy; 
the most significant feature is that when the emitting area takes the form of a line, 
it is shorter than for mercury at the same current. ‘lhe apparent width of the 
line for the alloy, as far as can be estimated, is between 10-* and 5 x 105*-ems 
leading to an apparent emission current density of 2:5 x 108 to 5 x 10° amp/cm’, 
whereas the apparent current density for mercury is only slightly in excess of 
10° amp/cm?. 

There is considerably less copious evolution of vapour from the alloy cathode 
than from the mercury cathode, and it may be that for mercury the large amount 
of vapour emitted renders the surface at the cathode spot indeterminate and 
lowers its effective work function. 

For the liquid alloy the emitting line is easily broken up into unit spots 
even when moving at a velocity of 104 cm/sec. For mercury the corresponding 
velocity is about 2 x 103 cm/sec. It is well known that the cathode spot exerts 
a considerable net downward pressure on the cathode under it, so that for the 
less dense liquid alloy rupture into dents (and hence unit spots) would take 
place more rapidly. 

All the photographs here reproduced are of a region of dense and highly 
excited vapour shot off from the cathode at the emitting area. It is not known 
how far out from the cathode this region extends, or if there is a ‘ dark space’ 
between this region and the cathode, as is the case for mercury (see Smith 1946). 
For mercury this ‘dark space’ extends outwards for about 10° cm., and as the 
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width of the emitting line appears also to be 10-? cm. for mercury, the writer has. 
suggested (Froome 1949) that the actual emitting line may be thinner than this 
and that a method of estimating its real size (as opposed to its apparent size) is 
to determine the maximum current which the spot will carry. For mercury 
the value is between 0-5 and 0-1 amp., and if at this value the spot is symmetrical, 
i.e. its length is equal to its breadth, the suggested real emission current density 
is 2x 10®to 107 amp/cm?. For the liquid alloy the minimum current lies between 
the same values; hence, by applying the same argument, the real emission current 
density may lie between the extraordinarily high values of 1:2 x10’ and 
6 x 107 amp/cm?, if we assume that for a current of 7 amp. the length of the 
emitting line / is given by /=Ki where K ~4-0 x 10-4 cm. amp. 

Even if we do not assume the real value of emission current density to be as 
high as that given above, and take only the apparent value of 2:5 to5 x 10°amp/cm?, 
we find that Langmuir’s (1923) original theory of field emission is well sub- 
stantiated without the assumption of an enhanced electric field at the cathode 
due to surface bumps. ‘This is the first time that emission densities as high as 
those predicted by Langmuir’s theory have been observed. 

Langmuir suggested that since electrons have a high mobility, the positive 
ion current to the cathode, across the cathode fall of potential, must be positive 
space-charge-limited. This means that a certain field must exist in order that 
this current can pass, and he assumed that this field would be high enough to 
extract the electron current by ‘cold’ or ‘high field’ emission from the cathode 
surface. 

Langmuir’s well known space-charge formula reads 


J 1 / 2e \12 V3? 
ile) E 


where J, is the positive ion current density, V the cathode fall of potential, 
d the thickness of the cathode dark space, M, mass of the positive ions, and 
e the electronic charge. 
To pass this current a field X is needed where 
M.\14 | 
X =47 (=) J NAV UA = (3-28 x 107)V2 44 volt/cm. ; 

M is now the atomic weight of the positive ions, and J, and V are in practical 
units. 

For the liquid alloy this reduces to X =2:6x104J,12, assuming for 
convenience a cathode fall of 10 volts. 

To obtain the electron emission density produced by this electric field, the 
equation of Fowler and Nordheim (1928) is best used. This states that the 
emission current density J_ for an electric field X is given by 

OZ 10 Xe ( — 6°8 x 10’d3? 


J 4 ) amp/cm? 


‘ p 
where 4, the work function, is 2 volts. 
Figure 12 shows the electron emission density plotted against various incoming 
positive ion densities. It is seen that as soon as the positive ion current density 
attains a value of 10° amp/cm? an electron emission density greater than this 
value is obtained. The results described show apparent emission densities 
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of up to 5x10® amp/cm?. At such a value a positive ion current density of 


10° amp/cm? is quite conceivable, so that the results described may be taken. 
to verify Langmuir’s theory. 
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Figure 12. 


Further, as for mercury, there is an apparent afterglow, lasting a microsecond: 
or so, in the wake of a moving cathode spot. 
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ABSTRACT. The field and lens constants are calculated for two typical magnetic lenses 
as used in electron microscopes, using the relaxation method developed by Southwell for the 
solution of potential problems. The results are compared with those based on simple 
analytic approximations to the field distribution in the work of Glaser and Ramberg, and 
are shown to agree closely as regards those characteristics on which the discussion of the 
performance of the lens is based. 


§1. INTRODUCTION 


lenses to have an accurate knowledge of the field distribution in the 

neighbourhood of the pole pieces. The determination of the characteristics 
of a lens, including the aberrations, can be carried out if the field distribution 
along the axis of symmetry of the lens is known. Bertram (1940, 1942) has 
solved the field equation analytically with boundary conditions approximating 
to those of a practical lens with unsaturated pole pieces, and finds a field expression 
-of the form 


|: is an advantage in many problems connected with the design of electron 


H =H, sech? az. 
‘Glaser (1941) has exhaustively studied the trajectories due to a field of the type 
jeden 
OTE GIae 


No attempt has, however, been made to solve theoretically the problem of the 
magnetic lens with saturated pole pieces, although this is the condition under 
which lenses are worked in practice. 

This paper describes the solution of the field equations by the relaxation 
method developed by Southwell and his collaborators (Southwell 1946) for the 
numerical solution of potential problems. 'The calculation has been carried 
out for both unsaturated and saturated pole pieces. 


§2, EQUATIONS FOR THE FLUX FUNCTION 

In the first part of this paper it will be assumed that the relation between 
the magnetic field B and the magnetizing force H is linear, that is, B=,H, 
where pz is the constant magnetic permeability. This is only true over a restricted 
range of H for materials used in practice, and the modifications required when p is 
not constant will be discussed below. 

The application of the relaxation method to a problem of field distribution 
is simplest when the field equations are expressed in a form containing a function % 


* This work formed part of a thesis submitted to the University of London for the Ph.D. degree. 
t Now at Royal Holloway College, University of London. 
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which is constant along the flux lines. is defined by the equations 


1 dus _ 1 op . 

> ror? HH, =~ a. als 08 8 @ (1) 
‘The equation for ¢ is : 

ay A 10/4 O0-4nnl 


See ed ek 2 ee (2) 


where mJ is the number of ampere turns in the coil and A is the total cross-sectional 
area of the windings. 

Physically, the value of 7 at any point gives the magnetic flux passing through 
a circle which lies in a plane perpendicular to the axis, with its centre on the axis 
and its circumference passing through the given point. 

Ataboundary between media with different values of » the tangential component 
of H is continuous, and the normal component of B is continuous. 

Thus at a point on the boundary 


(Hy), =(ies (HA (4ae nt (3) 
From equations (1) 


wizws (GS) =-CE) ee (4) 


ju On wu on 


\ 


§3. THE RELAXATION METHOD 


The solution of a potential distribution is given by the relaxation method 
in terms of a function defined by computed values, which satisfies the finite 
difference equation corresponding to the field at a finite number of points. The 
number may be increased indefinitely to obtain the degree of accuracy required, 
the only limitation being the labour involved in solving the system of simultaneous 
equations. The relaxation method is essentially a device for solving a large 
number of such equations quickly, and in such a way that mistakes can be detected 
and rectified at any stage. 


Figure 1. 


The points at which the values of the function are to be calculated are most 
conveniently taken at the corners of a square mesh, as shown in Figure 1;. The 
value of the function w, at any general point of the field may be found in terms 
of the values «,, ws, ws, w4, at surrounding points, and the mesh length h. Special 
equations are needed at boundary points between media of different permeability. 
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We require to solve the differential equation for the flux function over the: 
field of the lens. As the function will be computed in arbitrary units, equation (2), 
is replaced by 

C0 0a 1 Oe 


ve Oz? + or +7Z =(. 0, ©) 0 eens, (5) 

ys is then given by 
b  O-4nnl/A 
f = ——_ ere (6) 


where Z is a constant whose value depends on the units in which w is calculated. 

At ordinary points of the field the first-order finite difference equation 
corresponding to (5) is 

h h 
Wy, +(1 = z) Wo ty +(1 + 5) w,—40) +h?Z,r=0. 

This is obtained by expanding w in a Taylor series at the point O and neglecting: 
terms containing powers of h equal to or higher than the third. It is therefore- 
accurate to the square of the mesh length. 


The first step in the process of solution is to guess a distribution of w throughout 
the field. Then the residual 


h h 
Rymen+(1— 5) boy t(1 +) oy ten +8207 Be PD (7) 


is calculated for each point of the mesh. The residuals are now reduced to zero: 
by subtracting R/4 from w at the point where R is greatest, and modifying the 
- surrounding points accordingly. The process is repeated until all the residuals. 
lie between —2 and +2 (where the units of this residual are those of the last 
digit or decimal of w), after which no change in the values of w will reduce them 
further. Values of w which are fixed by the boundary conditions are of course 
never changed. Finally R is calculated from equation (7) for each point again, 
and any mistakes are detected by the occurrence of residuals greater than 2, and. 
can be rectified by reducing these as before. At the final stage the signs of the 
residuals should be distributed at random over the field, and the sum of the 
residuals should be almost zero. This reduces the possibility of cumulative: 
errors in any part of the field. The accuracy of the solution in any region can 
be increased by reducing the mesh length, or by taking account of higher order 
terms in the Taylor series when forming the difference equations. 

Special difference equations are required for points near a boundary where 
the distances between neighbouring points are unequal, the so-called ‘irregular’ 
stars’. Also it is convenient at some points to use a square net with sides. 
inclined at an angle of 45° to the z axis. This is required near diagonal boundaries. 
and in regions where a change-over to a net of smaller mesh takes place, and is. 
also useful for interpolating the value at the centre of acoarse mesh. At boundary 
points it is assumed as a first approximation that the value of 4 Was everywhere- 
of the order of 1,000, so that 1/ can be neglected in comparison with unity, 
giving a possible error of 0-1 °/, in the value of w. Then the boundary equation (4) 
reduces to (w|dn),;,=0. The region outside a sharp corner presents some: 
difficulty. The flux lines tend to concentrate round the corner, and the value 
of the flux function varies rapidly so that the Taylor series converges slowly.. 
If the corner were mathematically sharp the flux density would be infinite there,. 
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dut in fact the radius of curvature at the corner of the iron pole pieces is finite, 
and the mesh in this region should be at least as small as the radius. 
The difference equations for all these special cases are given in the Appendix. 


$4. CALCULATION OF THE LENS FIELDS AND TRAJECTORIES 

The flux distribution was investigated for a simple lens such as might be used 
for the objective of an electron microscope, the permeability being assumed 
infinite. ‘The boundary conditions of the problem are given by the results 
of the previous section, and the value of w along the axis is given by 


Os Ou 
= = —prH, = =yurH,. 
or _ dz at 
y% and w are therefore zero on the axis. 
At the first stage of the calculation the value of Z was not defined, but a 
reasonable distribution of w was assumed, and at each mesh point the function 


h h 
a +(1 — x) Ws toy +(1 + J) o1— hey 


or the equivalent expression for boundarv points, was calculated. ‘Then the 
average value of this residual was obtained and Z chosen so that the total residual 
over the area of the windings was zero. ‘This device was found to lighten the 
work considerably, since the only fixed values of w at which an accumulated 
residual could be ‘ thrown away’ lie on the axis at one corner of the field, opposite 
the air gap (marked X in Figure 2). 


=, 
= 


BES a ES 
==, 
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P Iron Casing and “ay 
Pole Pieces KN Winding Space 
Figure 2. Section across objective lens I showing flux lines in the windings. 


Rigorous justification of the method of solution is not necessary; the method 
is immaterial if the final values of w satisfy the original set of linear equations, 
which have a unique solution. The error involved in approximating to the 
differential equation by finite difference equations, and the possible accumulation 
of errors in the relaxation process itself, however, has to be discussed. It was 
mentioned above that the finite difference equations are an accurate representation 
of the differential equation if powers of the mesh length higher than the third 
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can be neglected. This can be tested in the case of any solution by halving the 
mesh length and finding the values of w at the new nodes. _If the values previously 
obtained in the larger mesh remain unaltered, the neglect of terms in the Taylor 
series is justified; if not, the mesh length must be further reduced until two 
successive calculations correspond. 

In unsaturated iron the permeability coefficient is about 1,000. Thus 
three-figure accuracy can be obtained in the values of w if infinite permeability 
is assumed. This accuracy can be maintained over the whole field irrespective 
of the absolute value of w, as may be shown by considering two sets of solutions, 


Uys tn, ay «atte 
and 
Uy +8Uy, Uy + 8g, Ug +8ulg, ... Uy +OU,, 


where 8u,,/u,<e and « is the possible proportional error of u, 1n one region of the 
field. The only fixed boundary values of w are those on the axis, where w is 
zero, and since the field represents a set of simultaneous linear equations in 
@ , Wy, ... W,, any set of solutions may be multiplied in the same ratio to give 
another set of solutions. ‘Thus 


Ou, _ lls éu, 


seee ) 


Uy Us te 


and the possible percentage error remains constant over the whole field. 

The absolute value of w falls rapidly towards zero near the axis, so that the 
values of w must be calculated to a greater number of decimal places near the 
axis than is necessary near the pole pieces (see Table 3), but since the possible 
percentage error remains constant, three-figure accuracy can be obtained through- 
out the field. Another possible source of error arises from the fact that the only 
fixed values of w lie on the axis, and errors in the rest of the field are likely to be 
cumulative. The sum of the residuals in any arbitrary area of the distribution 
should be numerically less than 10% of the number of mesh points in order to 
reduce this possibility. 

The absolute value of H, along the axis may be found in terms of the number 
of ampere turns in the coil by an application of Ampére’s law 


H.ds=0-47nl, 
re 
where C is a closed curve lying along the axis and round the outside of the iron 
case where the field is zero. Using equations (1) and expanding w in a Taylor 
series we obtain 


(HEELS) cacy oS CO) Re oe al ule (8) 
and 
Ae | Os 0 4als 
—~o Wm 
where «,, is the value of w(h) atr=h, z=0, and H,, is the value of H at r=0, z=0. 


The distribution of H, for lens I is shown in Figure 3. 

The calculation for this lens indicated that one of the flux lines (Y in Figure 2) 
is straight and parallel to the axis at the centre of the pole-piece gap. If this is 
assumed to be generally true for all lenses with similarly shaped pole pieces 
within the limits of accuracy required, the work of calculating the field distribution 
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is considerably reduced. The distributions of H, were found by this method 
for similar lenses in which the ratio of the gap to the inner diameter d of the pole 
pieces varies; these are shown in Figure 3, together with the distributions. 
assumed by Glaser (1941) and Ramberg (1942). 


% (cm.x 25) 


Figure 3. Variations of the field distribution with the pole piece ratio b/d. 


The equation for the trajectories of paraxial electrons in a magnetic lens is 
d*r —e 
eerste Pe ee 
dz? 8mc?V 0) 
A numerical step-by-step method of solution is suitable, as the values of H are 
available at intervals along the axis. Neglecting terms of order h*, equation (9): 


becomes 
é 


T,—2n +12 — 27,-n =f? Smc2V Vz} fabs. he 
w(h) eH 
i H,_,=H,— and ——, =k, 
Putting Z~h m iss Qmc2V 
,wr(h 
we obtain = {2 —kh? ( ’) hr | al 
Win, 


From this equation the radial distance of an electron is calculated at intervals 
z=h. w(h)/w,, is accurate to three figures, and five figures may be retained in 
the values of r,. Since the number of steps is never greater than fifty, the final 
values of r, may safely be regarded as accurate to three figures. 

To find the focal length, the initial path of the electron before it enters the 
field is assumed to be at a constant distance from the axis. ‘The focal length 
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is defined by analogy with geometrical optics. It is obtained by producing the 
final line of the path backwards to cut the initial path (Figure 4). _ sis the apparent 
shift of the centre of the lens system and f and s must be used in all calculations 
involving the formulae of geometrical optics. These formulae will only be 
valid if the object distance is greater than s. ; 

Figure 5 shows the focal length f plotted as a function of k=6-85n1/,/z, 
representing the power of the lens. 


S (cm) 
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| 6-85 
VY 
Figure 4. Figure 5. Variation of focal length with coil 
strength. 


§5. CALCULATIONS FOR SATURATED LENS 


The above results have been obtained by assuming that the permeability 
coefficient is greater than 1,000. This assumption is not valid for strong lenses, 
and its range of validity can be determined by finding the relation between the 
intensity of magnetization in any part of the iron and the number of ampere turns 
in the coil. In the lens under consideration it was found that saturation would 
occur if the number of ampere turns exceeded 1,200, and since the optimum 
strength of the lens has been found experimentally to lie between 4,000 and 
5,000 ampere turns, it is evident that under normal working conditions the core 
and pole pieces are likely to be saturated. The possibility of using the relaxation 
method to calculate the field of a saturated lens was therefore examined. 

If the iron in the field is saturated, the ratio B/H is not constant, and the 
relation B =H has to be replaced by 


B= f(H), 


where f is a function defined by the magnetization curve for the material. In 
an isotropic medium H is parallel to B, so that 


H,|B,=H,]B,. 


If is defined in terms of B as before, no simple equation can be found which is 
satisfied by #. It may, however, be assumed that the ratio B/H varies slowly 
over the field, and that in the region of any mesh point there is a linear relation 
between B and H. 'The coefficient j., is then defined as the slope of the tangent 
to the (B, H) curve. Thus B—B,=,,H, where By is the intercept of the tangent 
on the axis and ju, and By vary from point to point. 
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"The function % may be defined as above by the equations 


‘Then equation (2) follows as before, and from conditions (3) the boundary 


equations are 
Ow Ow 
@=@, and a bse an . 
where p= B/H. 


The process of solution is complicated by the fact that the values of u are not 
known until B has been calculated, and p» is involved in the calculation of B. 
A method of successive approximation is therefore necessary. Bearing in mind 
the physical conditions of the problem, a distribution of « along the boundary 
of the iron may be assumed and the flux function calculated as before, but using 
the exact expressions for the boundary conditions for finite u. Only the equations 
for the boundary points are affected by the variation of 4; those for the points 
inside the iron are exactly the same as for points in air (see Appendix). Given 
the flux function, B can be derived from equations (1), and hence new values of 
from the (B,) curve (Figure 6). After some experience of the behaviour of 
the distribution the computer can make the process converge fairly quickly. 
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Figure 6. Variation of magnetic permeability with induction. 


In order to investigate the field distribution when the pole pieces are saturated, 
the calculation was carried out for a lens similar to that illustrated by Hillier 
and Ramberg (1947, Figure 21) and used by them in attaining a resolution of 10 a. 
The distribution for unsaturated iron was first calculated. the field being obtained 
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inside the iron as well as in air. Table 1 shows the distribution in the coil 
windings and Table 2 that in the region of the pole pieces on a larger scale. The 
field inside the iron in Table 2 has been obtained by assuming a uniform distribution 
of flux lines across the sections marked AA’ and BB’ in Table 1, these sections. 
being taken at a sufficient distance from the pole pieces to avoid their disturbing 
effect on the uniformity. Table 3 shows the axial field distribution in more detail. 
The values of w in these figures are given in arbitrary units; the units in Table 3 
happen to be larger than those in Tables 1 and 2 by a factor 435/20. This 
calculation with infinite permeability showed that there is an area in the air gap 
where the flux lines are straight and parallel to the axis of the coil, that is, the 
line w=20 in Table 2. In such a region the radial boundary equation (see 
Appendix) with w,;=w, and Z,=0 reduces to 


h h 
(1 - ae (1 + J) 1 409=0, 


which is independent of ». ‘Therefore the introduction of a finite value of pw in 
this region does not affect the distribution, and it was assumed that saturation 
at the outer boundaries of the pole pieces could be regarded as having no effect 
on the region near the axis, and that the distribution of w along lines surrounding 
the saturated tips was unchanged. ‘Thus the problem was reduced to one with 
definite boundary values derived from the unsaturated lens. These boundary 
values are shown on the dotted line in Table 4. Nothing is assumed here about 
the absolute value of w along this boundary; this is undoubtedly changed by 
saturation at the outer edge of the iron; itis merely the constancy of the distribution 
along the boundary that is required. The absolute value is determined by Z, 
which is involved in the conversion of w into % (equation (6)). Since nothing 
is known about Z in this case, a different method has to be adopted in deriving the 
intensity from the function w, which means that there is an uncertainty about the 
coil strength for which the calculation is being done until it is completed. 

The process can best be explained by following the calculation step by step: 


1. B,, B, are defined by the relations 


K dw 
Beary 


where Kw=y%, and K in this calculation is chosen arbitrarily. The order of 
magnitude of K is obtained by considering the corresponding constant in the 
unsaturated case and increasing it proportionately to the number of ampere 
turns for which the calculation is required. 

2. A distribution of 4, along the boundary is then guessed, and the function w 
calculated over the field by the relaxation method, using the finite difference 
equations in the Appendix. 

3. dw/dx and dw/dr can then be calculated for points on the boundary by 
differences, and hence B, and B, from equations (10). Then B?=B,? +B2, and 
a second approximation to at each point can be read off the (B, ») curve. 

4. Steps (2) and (3) are repeated until B and uw correspond within the limits 
of accuracy required. 

The results for one coil strength are shown in Table 4. 
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The distribution of w along the axis is then known, and since K has been 
arbitrarily chosen, the distribution of y is also known, and the value of nJ for 
which the calculation has been done can be found from Ampere’s law as follows: 


. - (ew - 
From (8) PM GL. HITS Pea hs 2Ka(h). 
0 

r (2K ORs A 
Thus | He de=| R w(h) dz= ee 
where = (2 = 

Wm 
and w,, is the value of w(h) at z=0. 
But | H, dz =0-4nnl, 
so that nl = ne ¢ 

°L77 


When one calculation has been done for the highest value of J required, 
the other distributions can be obtained easily by choosing the initial approximations 
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Figure 7. Variation of focal length with coil strength and voltage. 


in the appropriate proportions. It was also sometimes useful to vary K in the 
middle of a calculation in order to reduce the residuals as rapidly as possible. 
The resulting coil strengths are shown in the table below. 


K 1-12 1-68 1:95 2-18 2°35 
On 1-045 1-14 1:20 1:29 1°31 
A 12-63 16-91 18-78 20-63 22°40 
Hy, (gauss) 5850 9570 11710 14090 15350 
nI (amp. turns) 1180 2580 3560 4620 5480 


When the coil strength was less than 1,000 ampere turns, it was found that 
the effects of saturation could be neglected, and, therefore, the distribution shown 
in Table 3 is valid for any coil strength below this limit. 

The focal lengths of the lens for various coil strengths and accelerating 
voltages were calculated by the method of $4, and Figure 7 shows the curves 
of focal length plotted against ampere turns. 


Table 4. Field distribution along the axis with a coil strength of 5,480 ampere turns. ‘The numerical values in italics are those of the total induction B, in kilogauss, the others are those of w. 
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As before, the accuracy of the relaxation calculation to the number of figures 
required was ensured by halving the mesh length until the distribution remained 
unaltered by further increase in the number of nodal points. This also justified 
the use of discrete values of the variable x at boundary points. The final accuracy 
of the results is limited by the knowledge available about the magnetic properties 
of the iron, but the value of ,« was not found to be very critical, especially in regions 
where w did not vary rapidly. A variation in » of the order of 10% usually 
only produced a variation of 1% in w. Thus if other sources of error are 
neglected, the error in w is not greater than 1%, and it is possible to find the 
focal length to an accuracy of 0:1 mm. 


§6. DISCUSSION OF RESULTS 


Comparison has been made above with the ‘classic’ field distributions of 
Glaser and Ramberg. Cosslett (1946) has based a discussion of the resolving 
power of the electron microscope on their work, and it is of interest to compare 
this with the results obtained by the relaxation method. 
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Figure 8. Variation of focal length with lens strength. 


The curves of f/d against a parameter representing coil strength have been 
derived for both these distributions and are reproduced in Figure 8. ‘The focal 
length fy in this case has been calculated according to a different definition from 
that given above. The alternative definition is given by Ramberg in terms of 
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a trajectory initially parallel to the axis, and the slope of the trajectory at the 
point where it cuts the axis. Thus fp =7/Ao, where Aj is the slope at the point 
r=0 (Figure 4). 

As can be seen from Figure 4, the two definitions are equivalent when the 
focal point is in field-free space. For strong fields fy decreases asymptotically 
to zero, while the geometrical focal length becomes infinite, indicating the 
existence of two or more focal points for fields of this value. Ramberg’s focal 
length is more convenient for the determination of the spherical aberration 
coefficient, although it assumes that the part of the field beyond the focal point 
does not contribute appreciably to the aberration. ‘This assumption is justified 
in practice, since lenses are not worked at strengths beyond the minimum 
geometrical focal length, and in the working region the two focal lengths do not 
differ greatly. 

In Figure 8 the abscissa kp? is the lens parameter as defined by Ramberg, 
and is given by 

hy? = PH,,2|(2-63)V. 

The Figure also shows the corresponding curves obtained by the relaxation 
method for the case of infinite permeability and for a coil strength of 5,480 ampere 
turns. ‘The agreement is close, especially for values of ky? between 60 and 100, 
and indicates that although the field distributions represented by Glaser’s and 
Ramberg’s expressions are not good approximations to the field of a practical 
lens, they do give a good idea of the variation of lens power with coil strength 
and accelerating voltage, and give a satisfactory basis for the calculation of the 
spherical aberration and the discussion of lens performance. 


S72 “CONCLUSIONS 

During the fifteen years or so of the development of electron microscopes, 
progress in their construction has been made chiefly by trial-and-error methods, 
guided by mathematical treatment where approximations made this possible. 
A great many empirical generalizations have been made as a result of this work, 
and optimum operating conditions are now almost standardized. In view of 
the constant testing that these generalizations undergo in practice, it seems 
unlikely that the application of relaxation methods to the standard designs would 
lead to anything new at this stage. The relaxation method might however 
provide a useful means of investigating new designs, and one which avoids the 
practical difficulties of direct measurement. 


APPENDIX 


The finite difference equations to be satisfied at various types of points are 
as follows: 
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Radial boundary points (Figure 9) 
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Cylindrical boundary points (Figure 10) 


hh h hh hp h 


Right-angled corner (Figure 11) 
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ABSTRACT. ‘The fluxmeter consists of a conducting system of small capacity which is 
alternately exposed to and screened from the electric field by the movement of a rapidly 
rotating metal shield, thus producing an alternating E.M.F. of period 0°83 millisecond. The 
-output, after amplification, is displayed on the screen of a cathode-ray oscillograph and the 
record carries an automatic indication of the sense of the field every 7:5 milliseconds. 

At the maximum sensitivity employed the device gives a deflection of one centimetre in a 
field of 20 volts/metre, with a background noise-level of 0:15 cm. (3 volts/metre). 

The response-time of the instrument is one half-cycle (0°42 millisecond). 

The fluxmeter has been developed to give the rapid response and high sensitivity 
required for studies of the electric fields of thunderstorms in the intervals between the 
-separate strokes of a lightning discharge. 


§L_ INTRODUCTION 

: HE electrostatic fluxmeter described in this paper has been developed to 
| provide a convenient method of measuring the electric fields produced by 
thunderstorms. The principle of the device was first described by 
Ross Gunn (1932) and it is sometimes referred to as a ‘field-mill’. It does not 
appear to have been developed before in a form which gives it either the 
‘sensitivity or the rapid time of response of the instrument to be described here. 
It has been used by us in studies of the electric fields of thunderstorms in the 
interval of the order of 0-03 second between separate component strokes, and it 
has an obvious application in providing continuous records of electric fields over 
long periods. With a time resolution of somewhat less than one millisecond, 
it is capable of resolving field changes which the beautiful capillary electrometer 
method of C. TT’. R. Wilson (1916) misses, and in our experience it is more suitable 
for field strength records of long duration than the combination of exposed 
‘conductor and cathode-ray oscillograph devised by Appleton, Watson-Watt and 
Herd (1926), particularly when studying nearby storms. The latter method 
-calls for high insulation resistance of the exposed conductor in order to maintain 
a long time-constant and this is a matter of some difficulty in the presence of 

heavy rain and point-discharge. 


§2. DESCRIPTION OF THE INSTRUMENT 

The instrument consists essentially of a conducting system of small capacity 
which is alternately and regularly exposed to, and screened from, the electric 
field by the movement of a rapidly rotating earthed metal shield. The 
exposure-screening cycle is carried out 1,200 times per second. 

The conditions under which the cycle is performed cause an amplifying system 
connected to the underside of the conductor to give rise to an alternating E.M.F. 
-of approximately sinusoidal form, whose amplitude is proportional to the strength 
of the field to be measured. ‘This £.M.r. is amplified and passed to the X plates 
-of a cathode-ray tube and the waveform is recorded by photographing the screen 
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in a drum-film camera rotating about a horizontal axis. The amplitude of the 
envelope of the crests and troughs on this record is proportional to the field 
strength and alters as the field strength alters. A special device described in §5 
indicates the sense of the field every ninth cycle and so shows which side of the 
envelope is to be selected for measurement. Some examples of records taken 
with the equipment using two different methods of recording are given in 
Figure 1 (see Plate). In all of these if the crests are enhanced in amplitude by 
“polarity pips’ the field is negative and vice versa. The arrows, therefore, 
show the side of the envelope which is used for measurement. 

Constructional details of the instrument are shown diagrammatically in 
| Figure 2, which is, however, not drawn to scale. 


Figure 2. 


The conducting system is formed by 18 metal studs mounted on a paxolin 
disc B. ‘These studs, four of which are shown in the diagram, are constructed 
from No. 4 B.a. cheese-head screws of diameter 6:5 mm. and are set in a circle 

of diameter 20 cm. ‘They are connected together below in the cylindrical brass 
box C and joined to the control grid of a low capacity triode valve V, (not shown) 
which is mounted inside C. The output of V, is led by cable to such further 
amplifiers as are necessary and from them to the cathode-ray tube. ‘The capacity } 
of the stud-grid system is 50 pr. 

The screening disc A is of brass 23-5 cm. in diameter and 1-3 mm. thick. It 
is provided with 18 equally spaced holes of diameter 2:0 cm. whose centres lie 
immediately above the studs. ‘The spacing between A and B can be adjusted by 
altering the position of B after loosening its clamping screws and is usually 2 mm. 
A is rotated at 4,000 r.p.m. by the steel shaft E, mounted on ball-bearings in the 
tubular housing D, and carrying a pulley which is belt-driven by a } H.P., 
1,400 r.p.m. a.c. motor. 

The box C is mounted, like B, on the bearing housing D by means of a collar 
and clamping screws and by moving it downwards one can get at the valve without 
dismantling the two discs. In order to provide the best possible shielding for 
this valve, the box is provided with an earthed cover in the form of a thin metal 
sheet attached to the top of disc B and suitably perforated to give insulation to 
the studs. 

The framework supporting the device is built of angle-iron with welded 
joints and carried on a horizontal pipe G, which enters a socket at the top of a 
vertical pipe H, one metre long, mounted on the roof parapet of the laboratory. 
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The whole framework is enclosed in a rectangular earthed metal case with the 
A side open and usually upwards. When it is necessary to run the instrument 
in rain or hail, it is inverted, so that A is downwards, and it is swung out over 
the outside wall of the laboratory so as to obtain a suitable ground clearance. 
Its sensitivity in the inverted position is reduced by a factor of two. 

An essential part of the equipment is the earthing brush F, which consists. 
of two springy phosphor-bronze strips pinching the shaft E between them. The 
surface of contact is slightly oiled and requires periodic cleaning. ‘The purpose 
of F is to get rid of frictional charges which collect in an irregular manner on the 
disc A as a result of its rotation and of the belt drive. The brush and the 
whole framework are connected to earth. Unless this earthing is very good the 
noise-level of these frictional charges swamps the output. 

To obtain a smooth output from the device, the discs A and B must be 
carefully constructed and quite flat. Any tendency to vibration in A must be 
eliminated since it gives rise to changes in capacity and so to an unwanted 
fluctuation in the output. 


§3F AMPEIEYINGSCiRCU LES 

Figure 3 shows the arrangement used in the head-amplifying circuit of V,. 
In order to obtain a low grid capacity, a Mullard E.F.37 pentode is connected 
up asa triode. ‘This valve was chosen for its low microphonic properties, and _ 
as an extra precaution it is suspended by rubber bands. ‘The effect of the grid 
resistor R, upon sensitivity is discussed in §7. The output from Vj, is taken 
from a cathode-follower circuit and is passed to the recording room through a 
screened cable 7 metres long. The H.T. and filament supplies for V, are also 
provided through screened cable from the recording room and the whole 
equipment can thus be operated from a distance. 


Figure 3. 


Further amplification as required is provided by the X plate amplifier of the 
Dumont 208B oscillograph used to display the final output. 


§4. METHODS OF RECORDING 
In the present work we have been concerned with the recording of field 
changes associated with lightning flashes and have been able to use special trigger 


methods and fairly high film speeds. ‘The drumcameraused by us hasa peripheral 
velocity of 44 cm/sec. 


An Electrostatic Fluxmeter of Short Response-time 405 


Two systems of recording have been employed. In the first, the time-base 
of the oscillograph is not used and the output is seen on the tube face as a 
horizontal line, which is screened off by a thin strip of black paper. An anticipatory 
lightning trigger (Schonland and Elder 1941) is arranged to operate in the initial 
stages of the stepped leader process of a flash to ground and to put a bias on the 
Y plates of the oscillograph before the leader has travelled very far. This bias 
causes the spot to move out from behind the screening strip and is kept on for one 
second so that the field changes taking place during the flash can be recorded. 
The method has the disadvantage that information about early first-leader 
field changes is lost, though it can be provided, if required, by another oscillograph 
kept continuously running. The first four records of Figure 1 are examples 
of this type of recording. 

In the second method no trigger is used but the adjustable time-base of the 
oscillograph is employed. This time-base is ‘locked’ to the alternating E.M.F. 
from the studs so that the unmasked screen would show a single stationary wave 
pattern whose axis is vertical and whose amplitude alters with the field. A 
paper mask with a narrow horizontal slit is attached to the face of the screen in 
such a way as to cut out all but the crest of this stationary wave. When the field 
is steady or is varying slowly, the crest can be kept hidden behind a small masking 
disc at the centre of the slit, by manual adjustment of the X shift of the oscillo- 
graph. When the field alters rapidly, the crest travels along the slit and its 
displacement is recorded on the drum camera. In the case of many lightning 
flashes the sudden changes of field due to the leader process initiate the movement 
required for recording. 

In Figure 1, SBC 4-1 and SBC 3-5 are examples of this type of recording; in 
the first the crest-spot has moved to a position above the masking disc, in the 
second it has moved below it. 


SS OAR TyYeeEN D1 GAWOR 
In order to show directly the sense or polarity of the electric field whose 
strength, as explained, is measured by the envelope of the amplitudes of the 
waveform, we have incorporated in the instrument a device which automatically 
increases the height of the crests of the waveform every ninth cycle if the field 
is negative in sense and the troughs every ninth cycle if it is positive. An 
inspection of the record immediately gives the required sense of the field. 


oes 
ORR: 


Figure 4. 


Figure 4, which is a plan view of a portion of the disc A of Figure 2, shows 
how this is done. HH are the exposure holes in the screening disc, with the 
studs of disc B momentarily below their centres. An extra pair of holes, KK, 
have been placed in the positions shown. Below each of these are two extra 
studs P and P’ fixed in slots in the paxolin disc B. ‘Thin wires join these four 
‘polarity studs’ to the output lead of the main studs. he phase of their 
contribution to the output can be adjusted by moving them in their slots. 
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When the disc A rotates in the direction of the arrow, the leading edges of 
the holes K and H uncover the peripheral studs and the studs P at the same 
moment, thus enhancing the output by the contributions from P. In the 
absence of studs P’, both crests and troughs will be larger than before. The 
P’ studs, however, become exposed when all other studs are being screened, so 
that their output is opposed in sense to that of the other studs. By adjustment 
of the position of the P’ studs their output can be arranged to neutralize that of 
the P studs during the half-cycle that they areexposed. In this way the ‘ exposure’ 
amplitude peaks of the output are larger than normal every half-turn of disc A 
but the ‘screening’ troughs are not affected. On reversal of the sense of the 
external field, the exposure half-cycle is that of the troughs and the polarity 
indication, as shown in the example marked SBM 1-4 (Figure 1) promptly shifts 
over to the other side of the zero line. All the records in Figure 1 show this 
polarity indication clearly. ‘The easiest to follow is SBM 1-4 which begins with 
a negative electric field as shown by the enhanced crests and then passes through 
zero to an almost equal positive value, as shown by the enhanced troughs. This 
reversal of field occupies 0-06 second. 

SBM 1-3, as shown by its enhanced crests at the start, begins with a large 
negative field which falls to nearly zero after five abrupt positive field changes 
caused by successive strokes of a lightning flash. In the record below the first 
the field is positive and five more strokes lead to a large positive field which ends 
the discharge. ‘The sense of the fields in the other three examples shown in 
Figure 1 is similarly identified. 

In practice it is necessary to adjust the height of the P studs with thin packing 
washers to bring them nearer to A and so to give them a slightly larger output 
than the P’ studs. 


§6. SENSITIVITY AND PERFORMANCE 

The instrument in the form at present used by us gives a peak amplitude 
of 1 cm. in a field of 20 volt/metre and thus easily registers the electric field 
of fine weather which is of the order of 100 volt/metre. We have been able to 
record with it field changes due to separate lightning strokes at a distance of 25 km. 
The noise-level on our present site is about 1-5 mm. at the maximum amplification 
used and at the worst times of day. In a quieter locality the sensitivity could be 
much increased by using a larger disc with larger holes and larger studs or by 
cutting sectors from the disc, as shown by the dotted lines in Figure 4, and 
mounting the studs on a circle near to the edge of the disc. Local sources of 
electrical interference have, however, prevented us from using a higher 
sensitivity. 

During rain when the arrangement is inverted and swung outwards, it gives 
half the above sensitivity when one metre from the building in a horizontal 
direction and 12 metres above the ground. 

The equipment is calibrated by taking readings in fine weather and comparing 
them with simultaneous measurements by the stretched wire method in a nearby 
open field: 

The electric fields of nearby thunderstorms are so large as to require reduction 
of the amplification provided in the oscillograph if the instrument is not to go 
‘off-scale’. This is easily met by adjustment of the calibrated attenuator in the 
oscillograph. It often happens, however, that there are charged clouds in the 
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neighbourhood when one wants to register small field changes due to distant 
storms and to use a high amplification. We have employed two methods to. 
balance out such disturbing fields. ‘The first makes use of a small insulated 
plate J (Figure 2), fixed to the box C and mounted above the rotating disc A. 
A controlled balancing potential can be applied to J to prevent the local field 
from driving the spot off-scale. ‘The second method is to insulate the whole 
framework at the socket G of Figure 2 so that it and the disc A can be charged 
to an appropriate balancing potential while the studs and cathode follower are 
separately connected to earth as before. 

Trouble was initially experienced with 50 c/s. hum from the a.c. mains but 
this was eliminated by balancing it against the output from a variable three-phase 
mains transformer. 


| $7. THE RESPONSE-TIME-OF THE INSTRUMENT 
Any sudden change AEF in the electric field strength must produce a change 
_—AQO in the charge on the upper surface of the system, given by 


—~AQ=kAAE/An, 


where k is a site-correcting factor and A is the effective surface area of the conductor. 
- Acorresponding opposite charge AO must appear on the under-surface. Provided 
the field change occurs during an exposure interval, the initial response of a device: 
of this kind is therefore instantaneous. It can be shown, however, that the 
final steady value of the amplitude is less than the initial response and that a 
certain interval of time is needed for it to reach this steady amplitude after a 
| sudden change. 

Let us assume for simplicity that both exposure and screening intervals. 
are equal to T and that both operations start and stop instantaneously. If 
Cis the capacity of the system, including the grid of V, and the connecting leads, 
and R, is the grid resistor of Figure 3, we define « as e-"/°™, It is then 
easily shown that successive half-cycles of exposure and screening give rise to 
changes in the charge on the under-surface of the conductor which are given by 
the terms of the converging series 


AO, —AOG(1—a), AQ(i —«+a*), ... AQ[1l—a+a?... +{—1)*ax"]. 


These terms converge to a final steady value AQ/(1+«) when z is large, 
so that the final steady amplitude is 1/(1+«) of that obtained from the first 
exposure ‘kick’, AQ. 

In our usual arrangement, Tis 0-40 millisecond, Cis 50 pr. and R, is 2 megohms. 
w is then 0-0183 and the final steady amplitude is 1-8°% less than the first ‘ kick’. 

The response-time of the instrument will be defined as the number of half- 
cycles, N, before the amplitude reaches 99%, of its final steady value. Since 
the 7th term in the above series can be written 


[1 +(—1)"a" Mt] AQ/(1 +) 


the number N is given by «***~0-01. 

For the usual arrangement employed by us we have «' =0-0183, «? =0-0033 
so that N should be one half-cycle. 

Though the actual conditions of operation depart from the simple assumptions. 
made in deriving this criterion, the instrument does, in fact, give a response-time 
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of one half-cycle with R,=2 megohms. With a period of alternation of 
-0-83 millisecond,and a response-time of less than one complete cycle, the device 
therefore satisfies the requirements for which it was designed. 

To obtain a higher sensitivity in observing distant thunderstorms we have 
-occasionally increased R, to 10 megohms. Since in this case e~°*/“* is 0-01, 
the response-time should be 8 half-cycles or 3-6 milliseconds. It is found 
in practice to be about 9 half-cycles. 

With this value of R,, « becomes 0-449 and the final steady amplitude is only 
0-69 of that given by the first half-cycle. 

The amplitudes of successive half-cycles, as can be seen from the series 
given above, are not at first equal on the two sides of the zero. ‘They only approach 
-equality after the lapse of the N half-cycles of the response-time. Thus records 
taken with R, =10 megohms, such as SAX2-6 of Figure 1, show an apparent 
shift of the zero when a sudden field change takes place. In practice this is some- 
times useful in picking out small changes of field. ‘The other records hardly 
-show the effect, since they were taken with R, =2 megohms and the zero shift 
after a field change lasts for 2 half-cycles at most. 
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ABSTRACT. A theoretical study is presented of the TM-waves that can travel along a 
slipping stream of electrons. A slipping stream is defined as one in which the electrons 
move in parallel paths with velocity which varies with distance transverse to the motion. 
It is found that amplifying waves can travel along a slipping stream for all frequencies. It 
is also found that the slipping-stream tube has the remarkable property of combining the 
characteristics of a two-beam tube and a travelling-wave tube. This is due to the occurrence 
in the stream under suitable conditions of resonance layers, which act as highly reactive 
impedance sheets and can guide TM-waves of slow phase velocity in the same way as the 
spiral or corrugated surface in a travelling-wave tube. 

The slipping-stream tube behaves differently according to the fractional velocity range 
of the electrons. If the velocity varies linearly from V_, to Vy across the stream and 
a=(V_g—V_ ,)/(Vq+V_,@)-<0-42 it behaves primarily as a two-beam tube with plasma 
resonance frequency w /2, and therefore having a cut-off frequency w,=wo/a/2 and 
maximum gain 2:1w»/V decibels per unit length, where wy is the plasma resonance frequency 
of the slipping stream; in addition it has a low gain for w>w,. However, if «>0-42 it 
behaves primarily like a travelling-wave tube and the rate of gain is about 0:53w,/ Vise 
decibels per unit length at all frequencies above the plasma resonance frequency. This 
tate of gain is low compared with the gain possible with a travelling-wave tube in which the 
electron velocity is V_, but it is achieved without the need for any external slow-wave 
waveguide. 

Wave propagation along a slipping stream inside a waveguide, which in the absence of 
the electrons can guide a TM-wave of slow phase velocity, is also discussed. It is shown 
that this slipping stream travelling-wave tube has very much the same characteristics as a 
travelling-wave tube with a uniform electron beam. Maximum gain occurs when the phase 
velocity vg of the wave in the empty guide is about equal to the velocity of the electrons 
nearest to the reactive impedance sheet of the guide—in the example considered this is the 
fastest electron velocity V,. A parametric set of curves is given relating the complex 
propagation constant to the frequency for different ratios vp/V 4. 


$1 INTRODUCTION 


N his original article on the electron-wave tube Haeff (1949) describes a tube 

in which a microwave signal can be amplified by interaction with a slipping 

stream of electrons. ‘The slipping stream is a cylindrical beam in which the 
electrons are constrained to move in an axial direction, by the application of a 
strong steady axial magnetic field, with velocity which increases with radial distance 
from the axis of the beam. The stratification of velocity is obtained by surrounding 
the beam with a cylindrical conducting sheath or drift tube. The electron 
space charge within the drift tube depresses the potential along the axis, so that 
outer electrons travel faster than electrons along the axis (Smith and Hartman 
1940). ‘This results in what we shall call a slipping stream of electrons. 

An analysis of the operation of this tube is not given by Haeff, although he 
does discuss a simplified theory of the two-beam tube in some detail. Experi- 
mental results are cited as proof that the slipping stream tube can amplify micro- 
wave signals. One of the most interesting characteristics of such a tube is the 
spectrum of frequencies for which amplification is possible, and this is not 
discussed by Haeff. 
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A steady slipping stream of electrons can also be maintained in an electron 
beam by the application of crossed electric and magnetic fields transverse to the 
axis of the beam. No strong axial focusing magnetic field is then required. 
Moreover when the electrons are perturbed by an oscillatory field they are no longer 
constrained to move in an axial direction but can oscillate transversely. However, 
to obtain parallel flow of the electrons in the non-oscillatory state the applied 
magnetic field must have a prescribed variation across the beam which depends on 
the electron density distribution, and this in turn is related to the velocity- 
distribution. If these values are not maintained the beam may alter in width 
along its length and the motion will cease to be single-stream. ‘The non-oscillatory 
state of this system is the same as the single stream or type-S state of a linear 
magnetron in the particular example in which the electron density and the 
magnetic field are constant. In the more general case it would correspond to a 
type-S state of a magnetron with non-uniform electron density’and magnetic 
field. Inthe oscillatory state, however, the periodicity of the field in the magnetron 
oscillator imposes the condition that the propagation constant must be real and 
the frequency, which is determined by the boundary conditions, must be complex 
for oscillations to build up in time, whereas in the slipping stream an oscillation 
of known real frequency impressed on the stream induces a wave to travel with 
increasing amplitude along the stream. The propagation constant, which is now 
determined by the boundary conditions, must therefore be complex. It is the 
slipping-stream tube with crossed electric and magnetic fields that we shall analyse 
in the following report. 

Bunemann (1944) was the first to givea theory of the build-up of small amplitude 
oscillations from the type-S state in a magnetron oscillator and to derive an 
instability criterion. ‘The analysis that we shall give for the slipping-stream 
amplifier will follow along much the same lines as that of Buneman; in particular 
we shall also use an action function to describe the electron flow. 

One of the main objects of our work is to determine the spectrum of the 
frequencies which can be amplified. We shall consider not only a non-enclosed, 
or free, beam but also one which is loaded along its length by impedance sheets, 
which in the absence of the beam can guide a wave of slow phase velocity. The 
latter represents the beam inside a slow-wave waveguide in the manner of a travel- 
ling-wave tube. Although a loaded guide is not necessary for supporting an 
amplifying wave in a slipping stream of electrons it is of interest from the point of 
view of coupling the beam to the impressed field at the input and of extracting 
the amplified wave from the beam at the output. Also if ohmic losses in the walls. 
of the guide are small it becomes apparent that a much higher rate of amplification 
is possible with a slipping-stream travelling-wave tube than with a free slipping 
stream. 


§2. SMALL AMPLITUDE THEORY 
(1) An Arbitrary Velocity Distribution 


Our object is to investigate the properties of wave motions of small amplitude 
which can travel along a parallel beam of electrons, wherein the mean velocity 
varies with position across the stream. We use a quasi-stationary approximation 
throughout, in which retardation effects are neglected.. For simplicity we 
consider only a two-dimensional variation in a cartesian coordinate system chosen 
so that there is no variation of any quantity in the direction of the x-axis. In the 
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non-oscillatory state the electrons are assumed to move parallel to the z-axis with 
a velocity which depends on y. ‘The stream is therefore stratified parallel to the 
z-axis. We shall further assume that in the non-oscillatory state the electrons. 
are constrained to move in parallel straight paths by the combined action of a 
magnetic field of vector potential (0, 4,, 0) and of an electric field EF. Then 
if V is the electron velocity, ¢ the electric potential, p the charge density, we can 
obtain a non-oscillatory state for which curl mV = eB by the use of the Hamilton— 
Jacobi equation (see Appendix I) and Poisson’s equation. For this state, 


V=[0, 0, V(y)], ] 
A=(0, 4,, 9], 
m 
$= 7 V0) =$o Pe | eee (1) 
_em[,aV (avy 
Piao dy a a) Wee 


Then an action function ys, of the electrons exists, and is defined by 
gradi) =mV — eA. 


Moreover, the condition curlmV=eB enables us to use the Hamilton—Jacobi 
equation, and is fulfilled if, as we assume in our case, the electrons leave the cathode 
with zero velocity and there is no component of magnetic field normal to the 
cathode surface (Appendix I). 


Therefore 
é ; 
oo = mV(y), | 
Cie ee) Bek be 525 stars (2) 
Bape tense. | 
oy ape y= 
and the magnetic field is [B,, 0, 0], where 
mdV 
= ee eS a ee eer 3 
a e dy ( ) 


The second equation of (2) defines the magnetic vector potential, from which 
the transverse magnetic field B, is obtained. 

It should be noted that both the charge—density distribution and the magnetic 
field, and hence the electric field E,, are determined by the velocity V(y). 

We consider a perturbation in which curl(mV-—eA) is still zero. ‘Then we 
specify a possible case if we satisfy 

gradys=mV—eA 
a ie align ey. Aes (4) 

and ep = — (6A + grad ys)? ye | 


(which ensure that the equation motion is ee together with Poisson’s 
equation and the equation of continuity. We consider a of the form 


b= ho +4n(ye", 


where 
6=wt —he. oe Cee ge (5) 
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In keeping with the restriction of the analysis to a quasi-stationary approxi- 
mation, any perturbation of the magnetic vector potential will be neglected. 
Inserting (5) into (4) and retaining only first order perturbation terms we get 


p=hothre™, 
1 hV 
where oy = S (1 = =) DSCC S gi Ta aera (6) 
; i Race Op. ORs sp 
Poisson’s equation is aye yas ps teste Mees ate shes 8S: (7) 


We assume a small perturbation re” on the charge density p. Then, from (7), 


a es -m(1 © =) Gar aA ees} | ee (8) 


‘The next step is to eliminate r from (8) using the equation of continuity. Thus 


J,=pV.=(p +16) (v- Tne") 


‘h | 
~pV + (Vr- “7 pt) gt ee (9) 
in (1 Ob 
and J, =pV y=(e +e") é te) 
Ob, , 
= EE arta (10) 
and the equation of continuity is 
es; 0d 0 Z 
aa By Reo Sole oc (11) 
Therfore from (9), (10) and (11) 
LV Ned : d[ dt, 
( a) == {-n pi, + = (oI. oe Wee (12) 


On eliminating + from (8) and (12) we get the following differential equation 
for py: 
E ( h “) d?V 
a fal Gl 5 
w w } dy* 
2 


ey 12 
a {(2) fe (1- yh Pa arihan htt,” (13) 


in which we have written 
em dy" 


(ii) Linear Velocity Variation 
We shall not make any attempt to find a general solution of (13) for %, but 
shall consider only a linear increase of velocity with y and put 
Vi = Vo Ov OOS ee ee eee (15) 
Then from (1) p=emQ?/e. ‘Therefore from (14) 
@) Wig! ot ie. over see ee (16) 


where wy is the plasma resonance frequency. 
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On substituting 


into (13), the equation for ¥, reduces to 
d*py 2x dy 
‘dx? * x#—1 dx 

Having found the equation that determines the motion of the electrons and 
the field inside the stream it is now necessary to ensure continuity of the field at 
the boundary of the stream. The next step is therefore to derive an expression 
for the radial admittance at the boundary of the perturbed stream. We use a 
method originally described by Hahn (1939). 

If the boundaries of the unperturbed stream are y = +a the boundaries of the 
perturbed stream oscillate about y= +a. For small perturbations it is sufficient 
to evaluate the magnetic field at the perturbed boundary y=a+A but to take 
this magnetic field as applying at the boundary y=a. In effect the perturbed 
boundary is replaced by the unperturbed boundary and the current between 
them is replaced by a surface current. 

From Maxwell’s equations we get an expression for H,. ‘Thus 


Cee rep ge es Me (18) 


oH OE, 
i nd Aa are ed terse Ny 

oz ot ‘ ae) 
Since H, =H, +e" ail CAD), 
and E,, = —0¢/dy, it follows from (19) that 

we Ody 2 
Se rhe POE Ae ah Ca tne tial eeah aeaet S ZA 
1 h dy ple ( ) 


Rinna 
where J, =,e". 


If A is the vertical displacement of the boundary we have, neglecting the 
displacement current, 


Hiefa-PA)—H fa\=JA=pV,A. wuss (22) 
Now, as shown in Appendix II, 
Se [iw pe ae to: (23) 


Hence on substituting from (21) and (23) into (22) we get 
wedd, 1. tpV Ax, | _hV) 


Fig +) — h oy hl mM ChY 


Pie Cones ee ey FAS (24) 
h dy mhi—(hV/w) oy - 
The second line of (24) follows from (10). 
Also E, = —0¢/0z =E,e”. ‘Therefore 
h hV 
Eanes (1 y ~) econ Gone ereer (25) 


Finally, the wave admittance in the positive direction of y is given by 


Sa ia iwe [wo\? 1 a, 
Mani ra ia o* (2) aaarjopay et 


Ww 
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If we normalize this admittance by substituting 


P=— POR ge (27) 
EW 
and change the variable to x using (17), (26) becomes 
ail L\. 1 dd; 
Pa)= 24 (1-3) o9 shellensy ss (28) 


We also require the value of the admittance of the field outside the stream at 
the stream boundaries. In the first instance we shall take the beam to be isolated 
in free space. ‘The field outside the beam must then fall off in amplitude with 
distance from the beam. It is sufficient to represent this field by a simple 
evanescent 'T’M-wave with phase velocity w/h in the z direction. ‘The admittance 


of such a wave is 
: € k 
vai) ey 


In the quasi-stationary approximation which we are using, h >k, so we take 


v=i,/(2\z. Tee: (29) 


Normalizing it as in (27) gives 
PQ) sl. We. aS eee (30) 


Similarly the admittance at y = —a is 
Resa) 2 ae Se Oe a ee ee (31) 


(11) Determination of the Propagation Constant for the Free Stream 


The propagation constant # is determined by matching the admittance P 
at the boundaries of the stream. 

Suppose we have found two independent solutions of the equation (18) for 
#4. Denote them by F(x) and G(x). Then the general expression for #, is 


(i=, CG) a ren ee dae eae (32) 
From (27) the admittance is 
1 1\ F’(x) +cG" (x) 
P a 
(x) tee (1 =) F(a) eClx nn, eke ee (33) 
and the boundary conditions are 
Dai) = — Pa = 1 a (34) 
where, by (15) and (17), 
One =2( v) 
Ca I Bivens 
9 M9 Vo (35) 
and rae -52(1- $F). ee 
m0 seed) Vo 


It will be seen from (35) that the unknown propagation constant and the frequency 
both appear in the expressions for the boundaries x,,4. ‘This complicates the 
procedure for evaluating h, 

Before we discuss the results of solving (33) and (34) we shall complete the 
analysis by considering a stream loaded with admittance sheets at the boundaries. 
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(iv) Linear Velocity Variation with Loaded Guide 


In order to study the interplay of the electron space charge and a slow 
"TM-waveguide wave we shall analyse the following system. ‘The mean velocity 


of electrons increases from 0 at y=—a to 2wya at y= +a. At the boundary 
y=—a we assume a perfectly conducting sheet and at y =a an impedance sheet 
with such properties as are necessary for the propagation of a wave of phase 
velocity vy between y = —a and y= +a in the absence of free electrons. 
For this electron-free wave 
E,=E sinh T(y +a), | 
1We 
H,=- =p Bo cosh T(y +4), pa ee er (36) 
¥(a)= ~— = coth (2aT), | 
where T? =h,?—k? and hy =w/vp. 
3 lwe 
For u)<c Y(a)=—- Ts COth (20n,) soap Seer es (37) 
0 
The normalized admittance at y=a is therefore 
P(x,)=- & COLA 205) eee nae ee Coes (38) 
0 
Consequently from (37), taking the velocity of electrons at y= —a to be zero so 
that Vy) =aoq, 
P(x_q) = 0 
h Zh 
and P(x,,)= — — coth (==), 
ho @o 


Let 7 be the ratio of the phase velocity of the wave in the empty guide to the mean 
velocity of the fastest electrons. ‘Then 


% $26 
=e Meee Bale ose Seen (39) 
VERILY Bien lh, Ve 
and from (37) h = hy _ =), Nhe (40) 
Ww 
where oe a 
og 


P(x_,) OOS 
ane Pe Rat eae (41) 
P(x,)=—n(1- 24 )eoth (=), 
1 1 \ F’(x) +cG'(x) 
“ OT es (1 % 3) F(x) +cG(x) 
From the first condition of (41) we deduce that 
pet aT a) 2a eee ee (42) 


G(x_a) : 
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Voth ARN Aa SOMES Dips) 

The analysis presented in §2 is based on a quasi-stationary approximation of 
which the main assumptions are (a) that the phase velocity of the propagated wave 
is much less than c, (5) that the effect on the electrons of their own magnetic field 
is negligible in comparison with the effect of the externally applied magnetic 
field, (c) that the oscillating field is of small amplitude compared with the p.c. 
field. In the results to be described in this section attention has been focused on 
only one mode of wave propagation. It is the lowest-order mode, which in most 
applications would be the dominant mode. There is, however, an infinite set 
of higher-order modes, some of which may travel against the stream and some 
may be amplifying. We shall not pursue their study further in this work. 


(1) Solutions of the Equation for the Action Function 


The key to the behaviour of the slipping stream as an amplifier is contained 
in the functions F(«) and G(x), two independent solutions of the wave equation 
for ,, equation (18). Since x depends on h, which is in general complex, we 
must study F and G for complex values of the argument. We have therefore 
tabulated F(x) and G(x) defined in the following way: 


if, Oya, BAT OF =0; 
G(x) =e F(x) logey(1—x) +9(x)], (1,0) =0, 
G(x) = F(—x), G' (x) = — F’(—x). 


Then F(x) =c,[G(x) loge,(1 +x) +9(—«)]. The property of these functions which: 
is of primary interest to us in this discussion is that F(x) and G(x) have logarithmic 
singularities at the points (—1,0) and (1, 0) respectively. If we refer to the 
definition of « at equation (17) we find that at the singularities, 


Wo 


= AV toe 


But the right-hand side of this equation is the plasma frequency w, increased by 
the Doppler factor due to its motion with velocity V in a medium in which waves 
travel with velocity w/h. The singularity therefore corresponds to a condition 
of resonance between the applied frequency w and the apparent plasma frequency 
measured by a stationary observer. 


Ww 


(11) Properties of the Lowest-order Modes in a Free Slipping Stream 


In the analysis and discussion of the results it is convenient to describe the 
velocity range of the electrons in terms of « =(V,—V_,)/2V>. When «=1 the 
velocity is zero at one edge of the beam and 2V, at the other edge. 

In Figure 1 we show the imaginary part of the propagation constant h as a 
function of frequency for three values of «. ‘These results were obtained by 
solving the equations (33) and (34) of §2(iui). The frequency axis is marked in 
units of the plasma resonance frequency wy and h is measured in units of Wel Vee 

Consider firstly the spectrum for « =0-25. Figure 1 shows apparently three 
amplifying and three attenuating waves, taking the minus and plus signs of 
Im (AV 9/9) respectively. (In the following discussion we shall confine our 
attention to the amplifying waves only.) The wave marked (1) has an amplifi- 
cation range 0<w<w,, and a dependence of gain on frequency which bears a 
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strong resemblance to that of a two-beam tube. Comparison of the two shows. 
that the gain—frequency characteristic of the slipping-stream tube is very nearly 
the same as that of a two-beam tube which has stream velocities V,, and V_, but 
plasma resonance frequency w,/2. Since the Larmor frequency w,, is equal to 
w,/2, it appears that the plasma resonance frequency of the two-beam tube, whose 
electron space-charge field is balanced by an electric field of positive charges, 
should be compared with the Larmor frequency of the slipping-stream tube, 
whose space-charge field is balanced by a transverse magnetic field. As in the 
two-beam tube the phase velocity is about equal to the mean velocity V of the: 


0:3 
=) A 1) 
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SJ 
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FF 
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Figure 1. Dependence of rate of gain on frequency in a free slipping stream. 


electrons. ‘Thus wave (1) amplifies by the same mechanism as the dominant: 
amplifying wave of the two-beam tube. It is therefore convenient to call it the 
two-beam tube wave. 

Waves (2) and (3), however, have no counterpart in the two-beam tube. On 
the contrary they are of the same type as the amplifying wave in a travelling-wave- 
tube. The effective reactive impedance sheet, which guides the slow wave, is. 
provided by a resonance layer within the stream itself. In order to appreciate 
this phenomenon the reader is referred to Figures 2 to 5. 

In Figure 2 the loci of the points x, and x_,, defined by (35), are plotted for 
a range of values of «. They are values of the variable x, which varies linearly 
across the beam from x_, at the low-velocity edge to x, at the high-velocity edge. 
In order to assist in the explanation by the following argument a sketch of the lock 
of x, and x_, for «0-25 is given in Figure 3. 
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Figure 2. The loci of the points xg and x_g for a range of values of a=(Vg—V_g)/2Vo. 
The numbers mark corresponding points of x,, and x-~g. 
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Figure 3. Sketch of the loci of x,, and x_, for a<0-42. 
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Figure 5 Sketch of the loci of Xq and x_g for 0:9 <a<1:0, 
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From the expression for x, namely 


it follows that the value of « at the point of intersection of the straight line joining 
xX, to x_, and the real axis is the frequency 022=w/w, and that 


Im (AV 9/wo) =3(Imx, +1mx_,). 


The straight line joining x, and x_, also defines the distribution of the action 
function ys, across the beam. Now, as indicated in §3(i) %, has logarithmic 
singularities at the points (—1,0) and (1,0). If, therefore, the line passes near 
to either of these points ,, and consequently £,,, will have a pronounced maximum 
at the level in the stream that corresponds to the nearest point of approach. 
Moreover it is known (Ramo 1939) from the study of wave propagation in a uniform 
electron stream of velocity V that a slow TM-wave can propagate along the stream 
provided Re(w—Vh)=+w,. In terms of x this condition is Rex=+1. It 
signifies the state of resonance between the applied frequency and the plasma 
resonance frequency as measured by a stationary observer. If the real parts of 
x_, and x, enclose the point x =1, or the point « = — 1, the resonance condition for 
the propagation of a slow 'TM-wave will be approximately satisfied at one level 
in the stream, which can be appropriately called the resonance level. ‘The 
resonance level behaves therefore like the corrugated surface or reactive impedance 
sheet of a travelling-wave tube: it can guide a slow wave. ‘The axial electric 
field falls off rapidly away from the resonance level in just the same way as from 
the effective impedance sheet of a travelling-wave tube. If therefore a resonance 
layer occurs in the beam of a slipping-stream tube we should expect the wave 
supported by it to interact with electrons in other levels in the beam in much the 
same way as the wave guided by the helix or corrugated structure reacts with the 
electrons in a travelling-wave tube. 

In the light of these remarks let us consider how the loci of «, and x_, move as 
the frequency is increased from zero. In Figure 3 x, and x_, lie at the origin 
when 2=@/w,=1. As Q is increased x, moves out along the curve OAB and 
x_, along OA’B’. For these waves the real parts of x, and x_, do not enclose 
either x =1 or x=—1 and so no resonance layer occurs in the beam. ‘This gives 
the two-beam tube wave of Figure 1. However, when x, reaches the point B 
its path forks. One branch moves along the real axis to C then up to the limiting 
point D, where x= —0-542 +0-0617, while x_, moves to the resonance point C’ 
and then out along C’D’. In the range BC the propagation constant is entirely 
real but it acquires an imaginary part beyond C, as shown for « =0-25 in Figure 1. 
‘The cut-off frequency at C is simply related to « by the formula 


Q, = (0-754 /.) +0:247. 


In the range C to D it is clear from Figure 3 that a resonance layer occurs 
in the stream. It appears at the low-velocity edge of the stream at frequency 0, 
and moves towards the high-velocity edge as the frequency is increased. ‘The 
phase velocity also increases, as shown by the curve CD in Figure 6, in keeping 
with the velocity of the electrons at the resonance layer. As {2 tends to infinity 
Im (AV,/w,) tends to 0:061/(1 +), which gives a gain of 0-53w,/V, decibels per 
unit length. 
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If we now return to the point B and follow the other branch we find that Xe 
moves along the real axis to the singularity at E and then along EF, while x_, 
moves to E’, where x =0-507, and up to F’, where x=0-542 40-0617. Figure 1 
shows that as «,, moves from B to E the frequency decreases to Q3 =(0-754/«) — 0-247 
and Imh=0, but as «, moves out along EF the frequency again increases and 
Im (AV 9/w) increases asymptotically to the value 0-061/(1—«) as Q> oo. This 
gives a gain of 0:53w,/V_,, decibels per unit length. It is evident from Figure 3 
that in the range EF a resonance layer occurs in the beam and that it moves from 
the high-velocity edge towards the low-velocity edge as the frequency increases 
from Q;. ‘The phase velocity also decreases, as shown by the curve EF in Figure 6, 


Limiting Value 
orfl=co 4 


W/W, 


Figure 6. The phase velocity as a function of frequency for a=0-25. 


We can now explain why wave (3) has a higher gain than wave (2) for Q >Q,. 
It is because the gain for a travelling-wave tube is larger the lower the phase 
velocity of the wave guided by its effective impedance sheet. In the slipping- 
stream tube this wave has phase velocity about equal to the velocity of the electrons 
at the resonance layer and for Q SQ, the resonance layer is near the low-velocity 
edge with wave (3) and near the high-velocity edge with wave (2). Therefore the 
gain of wave (3) should be higher than the gain of wave (2). 

The gains of waves (2) and (3) do not increase with frequency according to 
the law Im(hV /wo) «< w'® of the travelling-wave tube. That they should not 
do so is apparent when it is remembered that the slow wave guided by the resonance 
layer will gather energy only from electrons which have velocities slightly greater 
than the phase velocity of the wave. Slower electrons will absorb energy from the 
wave. In a travelling-wave tube all the electrons travel slightly faster than the 
wave. 

Having discussed the characteristics of the waves for «0-25 in some detail 
let us now discuss the effect of changing «. Referring to Figure 2 we see that as 
% is increased from zero to 0-42 the branch point B of Figure 3 moves from 
x=—0-715tox=—1. Whena>0-42 the locus of x, Splits into two quite separate 
parts as shown in Figure 4, Waves (1) and (3) merge into one without any 
overlap in the gain-frequency characteristic of Figure 1. The examples « =0:5 
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and «=0-8 of Figure 1 illustrate this effect. They also show the steep rise in 
gain at high frequencies as «> 1, according to the law Im (AV p/w) ~0-061/(1 —«). 
Wave (2) however, starts at Q,=(0-754/ax)—0-247 when x, is at E. Between 
Q, and Q, =(0-754/«) + 0-247 it has a real propagation constant and when Q>Q, 
it has a gain which tends to the limiting value Im(hV /w 9) ~0-061/(1 +«) as 
Q-+> o. 

That the behaviour of x, and x_, indicated in Figure 4 is not the whole story 
becomes evident when we consider what happens as « approaches the value 
unity. It is found that the loci of «, and x_, change in the manner shown in 
Figure 5. It now appears that the locus of x, has a third branch GHK PQ (for 
clarity we have omitted the branch ECD from Figure 5). Moreover there is a 
critical value of « for which branches OAF and GHKPQ meet. ‘This new branch 
passes through the singularity at x= +1 and runs along the real axis from x=1 
to x=0. When «a>0-946, Rex,>Rewx_, along most of the curve GHKP; 
then it gives reverse waves, of which only the attenuating ones given by the 
negative sign of Im/ are of practical importance. Amplifying reverse waves 
are also given by this theory but they must always be accompanied by attenuating 
waves of larger amplitude.* When «=0-946 branches OAF and GHK meet at 
the point S, where x,= —0-30 +1-987, and when «>0-946 they split into an 
upper branch FLG and a lower branch OMP. ‘The corresponding loci of x_, 
are indicated by dashed letters in Figure 5. 

The most interesting feature of Figure 5 is the continued rise of the curve 
LF as «—1. It can be shown that the curve LF tends to the asymptote 


Imx,=0-061(1 +a)/(L—«) as x_,>0-542+0-0612 


and this gives the high gain, Im (AV y/w,)) =0-061/(1 —«), as «1 at frequencies 
above w,. It should not surprise one that the gain tends to infinity on this theory 
as%—>1. As the limit is approached both the velocity of the electrons at one edge 
of the beam and the velocity of the slow wave guided by the resonance layer tend 
to zero. The system should therefore behave as a travelling-wave tube with 
vanishing electron velocity which, in theory, would also have infinite gain. 

It is worth while at this stage to compare the theoretical expressions for the 
maximum gain of the travelling-wave tube (T'WT), two-beam tube (TBT) and 
slipping-stream tube (SST). We obtain the following values, expressed in 
decibels per unit length: 


w\'3 w 
Gowr = (=) a ; 


Grpr = 43 a where V is the mean velocity of the two electron streams, 


Cae eel = for the TB'T-mode 


and 
Ga = 0:53 = 


—a 


for the TWT-mode. 


* A discussion of reverse waves will be found in a forthcoming paper by G. G. Macfarlane and 
Mrs. A. M. Woodward. 
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These expressions apply to ideal systems in which no metal and other losses occur 
and in which the electrons move in regular streams without any random thermal 
motions. 

Maximum amplification is obtained in all cases when the electron velocity is a 
minimum. It is therefore reasonable to compare the systems for the same 
minimum velocity. Moreover, in order to obtain gain at frequencies above wy, 


with the two-beam tube it is necessary for V_, to be nearly equalto V. Therefore 
for comparison of the performance of these ideal tubes as amplifiers for frequencies. 


much higher than wy we can put V= V=V_,. Then we get 


: = 6(w/w)¥8: 4-3 : 2-1: 0-53, 


2 


A more detailed discussion of the effects of losses, thermal motion of the electrons, 
and possible feeding systems is beyond the scope of this paper. 

It still remains to mention another result which may be important in consider- 
ation of the use to which a slipping-stream tube might be put. In Figure 5 it 
is shown that the portion of the real axis from x = 1 to x = «© is acommon locus for 
x, and x_, for all values of x. There is no gain or attenuation for these waves and 
they are reverse uniform waves. When « is equal to unity the phase velocity 
depends on frequency for Q>1 as shown in Figure 7. These reverse waves may ~ 


OQ» 10,2 -20h © 30) £40 S060 G70 ae coe en Soleenico 
W/W, 


Figure 7. The phase velocity of the reverse uniform waves as a function of frequency, a=1. 


constitute a serious drawback to the use of a slipping-stream tube as an amplifier 
unless it can be found how to avoid exciting them. The fact, however, that Haeff 
has obtained a stable amplifier using a slipping stream is encouraging. 
Summing up the discussion we can say that a slipping-stream tube will behave 
differently according to the fractional velocity range of the electrons: if 


Cee an Vaie VW, =e VS) =a4<0-42 


it behaves primarily like a two-beam tube with cut-off frequency Q, =1/a./2 
but it has only about one half its maximum gain, viz. 2-1w/ V, decibels per unit 
length. In addition it has a gain of 0-53 /V_, for Q >Q,. However, if 
%>(-42 it behaves primarily like a two-wave tube with a gain of 0:53w/V__ 
decibels per unit length. Comparison with the travelling-wave tube and two-beam 
tube on the basis of the same minimum electron velocity shows that the slipping- 
stream tube has a lower maximum rate of gain than either of the other two. 
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(iit) Properties of the Amplifying Wave in a Slipping-stream Travelling-wave Tube 


In this section we give results for a slipping stream inside a waveguide. ‘The 
waveguide is such that a slow TM-wave of phase velocity vp can travel along it 
when there is no electron beam inside it. The velocity is taken to vary from 
0 to 2V)(«=1). The system is a form of travelling-wave tube using a slipping 
stream of electrons instead of a uniform stream. 

The theory, by means of which the following results were obtained, is set out 
in§2(iv). It was shown in that section that the waveguide could be defined in terms 
of an impedance sheet. This concept obviates the necessity of specifying the 
physical form of loaded guide required to support a slow wave of ‘M-type. 
It effects at once a unification and simplification because a TM-wave of phase 
velocity vy with the field distribution given in §2(iv) can be guided by different 
physical structures. The object of our analysis is to discover how the propagation 
constant / varies with frequency for a range of values of vy in the neighbourhood 
of 2V,. In most physical structures used to obtain such slow waves the phase 
velocity v, will vary with frequency. By giving a map of h marked in contours of 
constant 2V4/v,» and w/w, the dependence of the propagation constant on frequency 
can be quickly determined when the variation of vy with frequency is known. 
Such a contour map is shown in Figure 8, in which, instead of h, the quantity 


1-4 2% 


Figure 8. Contour map of the propagation constant for the slipping stream travelling-wave 
g tube vy is the phase velocity of the "[M-wave in the empty guide. 

given is x,=(w/w)—(2AV /w9). Since «=1, x_,=w/w9, so that the real 
component of / can be calculated immediately from the curves. The ordinate 
is simply Im(2AV,/w,). From the figure it is therefore clear that the maximum 
rate of gain is obtained when vp=2V 4, that is when the phase velocity in the empty 
guide is about equal to the velocity of the fastest electrons. It is interesting to: 
compare these results with the corresponding results for a travelling-wave tube 
with uniform beam, which are shown in Figure 9 (Woodward, unpublished work). 

In this figure the velocity of the electrons is V and the phase velocity in the 
empty guide isv. ‘The similarity of the characteristics at any rate for w>wy is. 
immediately evident. ‘The maximum gain is slightly greater for the slipping-stream 
tube. ‘There is, however, one noticeable difference. When v)<V there is an 
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upper limit to the gain with the travelling-wave tube, but when v»><2V, there is 
no upper limit with the loaded slipping-stream tube. The reason for this is that 
the electron velocity ranges in value from 0 to 2V, in the slipping-stream tube so 
that there are always some electrons for which the condition, v)>V, for gain at 
the highest frequencies is satisfied. However, for vy>2V, the gain falls to low 
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Figure 9. Contour map of the propagation constant for a uniform stream 
travelling-wave tube. 


values at high frequencies because the electrons which can give up energy to the 
waves are further from the impedance sheet. Moreover the field falls off rapidly 
from the impedance sheet. 

In Figure 10 we have plotted the phase velocity of the wave as a function of 
frequency. It is always less than the velocity of the fastest electrons. 


Figure 10. ‘The phase velocity v as a function of frequency for the slipping stream travelling- 
wave tube. The parameter is 2V,/v), where vy is the phase velocity in the empty guide. 
§4~ ADDITIONAL NOTES 
(a) ‘The non-oscillatory slipping stream requires applied electric and magnetic 
fields to ensure parallel motion of the electrons. From (3) the magnetic flux 
density is vs 
fog = 5 o> 5:68 x 10-8w» gauss. 


From (1) the potential in the beamis¢=mV2/2e. If the electric field is maintained 
by a potential difference between two parallel plates at y= +6 the potential on 
the plates must be 


m m 
aes 2e ee) -: V2; $(—6)= = Ve—(b-a)“y,, 


é é 


where ¢=() is the potential of the cathode from which the electrons are emitted 
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(6) In order to study the basic phenomena we have described a planar form 
of slipping-stream tube. Other forms can be thought of which may be of greater 
practicability. For example an annular beam might be used between the inner 
and outer conductors of a coaxial line; the magnetic field, which would be 
azimuthal, could be produced by axial currents in the inner and outer conductors. 

(c) From the small-signal theory it is not possible to deduce the behaviour of 
the system at large signal amplitudes. Nevertheless it is interesting to try to 
discuss a few general points connected with the transformation of energy in the 
system. 

It can be easily shown that the mean flow of energy across the section of the 
tube is 


— 


Wee = ib (Ex). nda- ~ { (VAG ndgne ee (43) 


This latter integral is the mean rate at which kinetic energy of the electrons is 
passing the section [. 

At the input section the R.F. energy can be neglected, V=V)+woy, and 
T extends from y= —a to y= +a. Evaluating the integral we get 


Ww= = pa( Wee + aw") Ve Gutheieh B10 (44) 


At a section farther along the tube this power is partly converted into oscillatory 
power. If we still think of the beam in two parts, one D.c. and the other a.c., we 
see that the A.c. power can increase at the expense of the D.c. power either by a 
decrease in mean velocity or a decrease of density of the electrons or by both. 
Now in a travelling-wave tube with strong axial magnetic field the D.c. power 
comes primarily from a decrease in mean velocity because lateral spreading is 
inhibited, but in the slipping-stream tube the beam is at liberty to spread and lose 
average density as well as mean velocity. On this account it seems likely that there 
will be a more efficient transfer of power from D.c. to A.c. at large amplitude with 
the slipping-stream tube. 
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AUP PEND I Xeel 


DERIVATION OF THE HAMILTON-JACOBI EQUATION 
The following discussion is restricted to non-relativistic motion. 
In terms of the scalar potential ¢ and the vector potential A the equation of 
motion of an electron in Newtonian form is 


mT =~ E-eVB=e(F +-gradp—Vacurl A). “ogee (45) 


For electrons moving in a stream the left hand side of (45) can be transformed, 


as in classical hydrodynamics, to give 


dV ov ON cae 5 
age mae +(V.grad)V = ai +4 grad V?—VacurlV. ...... (46) 


Then (45) becomes 


2 (mV —eA)=Vacurl (mV —eA)+grad(e6—4mV?). ..... (47), 


If the motion is such that 
Va. curl (nV. = cA) = 05 eee (48) 
mV —eA must be the gradient of some scalar, % say. ‘Then 


mV — 6A= gradds > a0". = Se > yates (49) 
and, on integrating (47), 


ee +imV?—eb=0, 7 ae See (50) 


since the constant of integration, which may be a function of time, can always be 
reduced to zero by suitable choice of ¢. 
Eliminating V from (49) and (50) gives 


S + 3m(eA + grad #)? =ed, 
which is the Hamilton-Jacobi equation in §2 (i). 

It can be shown (Gabor 1945) that the condition curl mV =eB is satisfied at all 
points of the motion if the electrons leave the cathode with zero velocity and there 
is no component of magnetic field normal to the cathode surface. 

Given that curlmV=eB, the two-dimensional problem treated in this paper 
demands that the magnetic field maintaining the velocity distributions be at 
right angles to the plane of the motion, and, in addition, be of a specified magnitude. 

For if V=[0, V,,, V,] and V is independent of x, then 


curl V= Ee - Sos 


so that B=[B,, 0, 0], where 
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APPENDIX II 
If the equation of the boundary surface is 
F(y, 3, t)=y —a—Aei—™ — 0, 
the function F must satisfy the equation (Ramsey 1935) 
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Theory of the Production of an Ionized Layer in a Non- 
Isothermal Atmosphere Neglecting the Earth’s Curvature, 
and its Application to Experimental Results 
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Rhodes University College, Grahamstown, South Africa 


Communicated by M. Blackman; MS. received 31st March 1949, and in amended form 
26th Fanuary 1950 


ABSTRACT. A new theory of ionospheric layer formation is developed, in which the 
temperature is assumed to vary linearly with height. "The equations are compared at each 
step with those obtained by Chapman in his theory of layer formation in an isothermal 
atmosphere. 'The equations for the maximum of electron density and its height are also 
given. The effect of the parameters on the shape of the layer 1s shown in graphical form. 

The equations are somewhat complex in form, but an ingenious graphical method has 
been devised suitable for the application of the theory to results given in the form of a 
Booker and Seaton parabolic distribution of electronic density with height. By application 
of the theory to mean hourly values for four summer months in South Africa, values are 
obtained for the temperature gradient, the temperature at 200 km., and its variation over the 
middle part of the day. The results obtained are in accordance with previous estimates and 
offer numerical confirmation of the theory that the atmosphere expands bodily upwards 
during the middle part of a summer day. 


Sail, ANS, GSU aON we 
1 (i). Introduction 


H£ most prominent among the theories put forward to account for the 
formation of ionospheric layers is that of Chapman (1931, 1939), which 
attributes their production to the ionizing action of solar ultra-violet 

radiation on the earth’s atmosphere. Other possible causes, such as corpuscular 

radiation, cosmic rays, alpha and beta rays from various sources, meteors, etc.,. 

are entirely neglected; owing to the irregular nature of most of them, thev are 
2 G-2 
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intractable to theoretical treatment at present. In Chapman’s theory a large 
number of simplifying assumptions is necessarily made: 
(i) the radiation is monochromatic; 
(11) che density of the atmosphere, on which the radiation is incident, varies 
exponentially with height ; 
(11) the atmosphere is at rest and uniform in composition ; 
(iv) the temperature is the same at all points in the atmosphere at all times; 
(v) the force of gravitation is constant with height; 
(vi) all the absorbed energy causes ionization. 


‘Charts showing the ionization of the upper atmosphere have been constructed 
by Millington (1932, 1935), using Chapman’s theory, but do not seem to have 
had much application. 

Various approximate theories have been put forward by Hulburt (1928, 1938); 
finally an attempt was made (Hulburt 1939a) to fit a theoretical curve to the 
observed mean values of maximum electron density for the E layer. 

Comparison shows that Chapman’s and Hulburt’s theories are the same in 
fundamental content, being based on the same initial assumptions. Their 
application to the E region has met with considerable success, but few attempts 
have been made to apply them to the higher more irregular regions. 

To develop a theory which will represent more nearly existing conditions - 
in the F regions, it seems advisable to reject Chapman’s assumption of an 
isothermal atmosphere and, instead, to allow the temperature to vary with height 
in a known way; as the simplest variation possible, it has been assumed in this 
paper that the temperature varies linearly with height from a given level upwards. 
‘The distribution of molecules with height in such an atmosphere will no longer 
be exponential and will first be derived. 


1(i1). Distribution of Molecules with Height 


Assuming no vertical movement of the atmosphere, the increase in pressure dp 
corresponding to an increase in height dh is given by 


— dp = gp dh, 


where g is the acceleration due to gravity and p the density at height k. We use 
the simple gas laws, and put v’ for the number of molecules per cm? at height h 
and ny for the corresponding value at a reference level o, the temperature at this 
height being 7). On substitution and integration between the limits mand n’, 
corresponding to heights 4, and h, we have 


; “h 
n' =n) exp ( = oe | 7) EL acs & (2) 
~ ho : 


where k= Boltzmann’s constant, and m=mean molecular mass. 


In this paper we limit the discussion to the case where TJ is a linear function 
of height and is given by 


f=T,+(h=}): —— 2 Se (3) 


Substitution in equation (2) and integration then yields 


; ) h —hp —(1+C/y) 
n=] by ’ ona ely 


0 
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where we have written 


ye 
oh ee ae Lave (5) 


On expressing this as a binomial expansion and inserting the limit y=0, 
equation (4) reduces to 


, ue 
n =niexp | - Fh}, Pen ree (6) 
0 
which is the exponential distribution assumed by Chapman and Hulburt. 


Table 1. Values of n’/nj for various y’s 


h n’|no 

(km.) y=0-00 ye) 2 405 y=8:10 
(deg/km.) (deg/km.) (deg/km.) (deg/km.) 

200 1-00 1-00 1-00 1-00 

225 3-64 x 10+ 3-44 «10-4 3°25e< 10m DDD NO 

250 1ossect Om 1i-32p<10q2 (SOO ee 2a el Oat 

275 4:78 x10? 5:54 10-7 5-98 x 10 62 52610—2 

300 1:74 10 25S el Oe S-O55-4\ Ome 3:64 x 107 

325 SQ SO 22 Ome 1695110 =e DOMES MOS? 

350 DADS AM 6:24 x 10-3 QB KO oS S< O—= 


The effect of a temperature gradient y on the molecular density may 
be seen from Table 1, which has been worked out assuming the values 
h,=200 km., 7)=400°K., g=900 cm.sec?®, m=15x1-66x10-4 gm. and 
m1 -38x<10-* erg. deo. 


1 (iii). Intensity of Incident Radiation 
The intensity J of the incident radiation (assumed monochromatic) will 
decrease as it passes through the atmosphere. Suppose that the intensity is. 
increased by an amount d/ in increasing the height from h to h+dh. If x, the 
solar zenith distance, is small, so that we may consider the earth’s surface to be 
flat, 
Glan see. die mr 2 6 ee es (7) 


where f is the atomic absorption coefficient. Substitution from equation (2), 
and subsequent integration yields 


rh { h dh 
lat, exp {mj TyB see x | pee (- ey 7) an eee (8) 


in which J,, is the intensity of the beam before it has reached the earth’s atmosphere 
and / the intensity at height A. 


1 (iv). Rate of Ion Production 


Suppose g is the number of ion pairs produced per cm? per second at height h. 
Then if w is the energy absorbed in ionizing one molecule, and assuming that 
all the energy absorbed produces ionization, 


Fee ie a a Se ge (9) 


W 
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Hence, on substituting from (8), 


h i h dh 
= file 2, exp{ - aie dh as Tepsec | pon(- Cc Tyan. 


With a linear temperature gradient, according to equation (3), this becomes 


ha=hA\=ee 
oe secx (1 +y 7) f 


I= 4%) 


a ee 
Saas ey ape exp 


oy aay (11) 
On proceeding to the limit y =0, as before, this reduces to 
W=% 
j € 
= Pe exp | — 5 (h— hy) — Mir sec xexp| — 7 (hh) | eee (12) 
0 


which, in the present notation, agrees with Hulburt and with Chapman. 


1(v). Layer Formation 


If ions are produced according to equation (11) or (12), there will be an equal 
number m of positive ions and electrons per cm3. ‘The rate of increase of ion 
density is then given by the well-known equation 


dn 


oF al eal ene Se, (13) 


in which « represents the effective recombination coefficient. 
Case (a).—It is instructive to consider first the solution of equation (13) for 
the constant temperature case. Substituting from equation (12), 
dn as ws ni, ( nBT, G H} 
— - — — (h—h,) —-—— ——(h-h ; 
Paar te, ex Pp T, (h—hg) G Secx exp T, ( 0) i 


According to Chapman, this equation cannot be solved explicitly. Hulburt, 
however, has pointed out that during the middle part of the day dn/dt is small, 
and may be assumed to be zero to obtain an approximate equation. ‘This 
assumption is confirmed by the experimental results to be described in § 2 of this 
paper. Equation (14) then reduces to 


(O 1 G 
n.2 = -_-—_ — — — =U) leg, oe neta eaeree 1 
nN, Bexp{ T, (h—hy) — FexP | T, (h hy) | \ ae (15) 
where we have written 


B= Png l slaw ve Wn Satie ard eee (16) 


and 


F=C cos x/nBT>. vie melee (a 
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Case (b).—A linear temperature gradient, given by equation (3), will now be 
considered. As the form of the equation is more complex than in case (a), we 
proceed at once to make the approximation dn/dt =0, obtaining 


A h—h,.\-@+¢”) 1 h—-h\-c 
m= B(14y 7) exp | — p (147 7) \. ahr (18) 


Equations (15) and (18) give the value of the electron density at any point in 
the layer. 


1 (vi). Maximum of Electron Density 


The maximum-with-height value N of the electron density, x, may be found 
by equating dn/dh to zero. When this is done with equations (15) and (18) in 
turn, we find that the required conditions are respectively 


1 G 
pep, =F (th) } te es eee eee (19) 
and 
1 h —hy =O y 
a(lty T, ) =l+G- epee (20) 
When substituted back in the original equations, we have 
DID Tex) (a= 1) ean Me ae pee (21) 
and 


N2=B {F (1 zs x) im exp { x (1 a 2) } Bt een (22) 


giving the required maxima for the two cases. 
The height of maximum electron density, /y, is given by equations (19) and (20), 
which, for future reference, are now written as 


(Ay), =o — = iy eae Se eee ie ever (23) 


(iy =ho + 2] (1 rZyrn -1]. Shae (24) 


1 (vii). Comparison of the Two Cases 


The changes in behaviour brought about by the assumption of a linear 
temperature gradient with height may be illustrated by means of the ratios 


Np (hy) —Ng 
—— and Se EO 
Na (h Ma —hy 


By this method the use of a number of questionable constants may be eliminated. 
The theory assumes that the composition of the atmosphere at the reference 
height hy should be the same as that throughout the layer. ‘There is reason to 
believe that the atmospheric gases exist mainly in atomic form from 200 km. 
upwards (Hulburt 1939b); Ao has therefore been chosen as 200 km., because the 
value of 7, becomes less certain as the height increases. Further, putting 
m= 15 atomic weight units, g =900 cm.sec® and k=1-38 x 10-6 erg.deg™, the 
value of C becomes 16:2 deg.km-!. Finally, for purposes of this comparison 
only, assuming 7) =360°K., m= 10" cm-$ and B=10-"’ cm’, on very debatable 
grounds, and putting y=0 for simplicity, we find F =4-4 x 10™. 
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Dividing equation (22) by (21), we obtain 


N,\2 air y (1+y/C) _yIC 25 
(3) = FY (14h exp (y/C)S eee (25) 


and from equations (23) and (24) 


(1 “ 1 

+ at AEs p= 

Ure hates CS es C7 ee (26) 
(Ay) Ko MA In F 


To obtain some idea of the absolute difference in height of the maxima, we use 


Una ee Ei Et ‘: 4 ies ee i} “ 2, In FI, Cer (27) 


obtained by subtracting equation (23) from equation (24). 
The values obtained are shown graphically in Figures 1 and 2. These graphs 
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Figure 1. Figure 2. 


show that a positive value of y produces a lowering of the maximum electron 
density in the layer and raises its height. 


I (vin). Effect of Variables on the Shape of the Layer 


As the general theory, summarized in equations (18), (22) and (24), does not 
lend itself to easy physical interpretation, it has been thought advisable to represent 
graphically the variations in shape of the layers brought about by changes in the 
variables. In order to eliminate as far as possible numerical errors arising from 


the use of uncertain values of the constants, equations (18) and (22) may be 
divided, giving 


my \2 sees y\ ate) 
) =(14y Th iF ms: 


enn —Cly 
xexp{(142)-5(1 typ) \ htaten (28) 


in which B does not appear. On the right-hand side of equation (28) we are 
thus left with four parameters, the height f,, the temperature JT, at 200 km., 
the temperature gradient y, and F. In Figures 3, 4 and 5, »,/N, is plotted against 
h, for various values of y, Ty and F, using equation (28). 
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Effect of y on shape of layer. 
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Figure 5. Effect of F on shape of layer. 
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In Figure 3, the curve for y=0 is drawn according to Chapman’s theory, 
and represents the fractional distribution of electrons at constant temperature. 
We note that a positive y raises the height hy. In addition, it is seen t ) increase 
the thickness of the layer, and the electron density above the maximum decreases 
more slowly than in an isothermal atmosphere. The importance of a knowledge 
of true layer thickness in any attempt to estimate y is thus clearly brought out. 

From Figure 4, it will be seen that an increase in T, produces an increase in 
the height of maximum electron density, but the thickness of the layer is not altered 
to any great extent. In Figure 5, the thickness of the layer is again practically 
unchanged by a change in F, but the height 4, decreases with increase in F. 
Remembering that for a given Ty (and therefore nj) F is proportional to cos y 
from equation (17), these curves show a result analogous to the theory for an 
isothermal atmosphere, i.e. the height of maximum electron density decreases to 
a minimum at local noon, all the other parameters remaining constant. However, 
since the theory of expansion demands atmospheric heating during the middle 
part of the day, Ty and y cannot be taken as independent of time. If the increase 
in Ty is greater than that of cos x, it alone would account for a rise in hy, even 
if y remained unchanged. 


$2. EXPERIMENTAL DATA 


Before developing a method by which the theory just expounded may be 
applied to experimental results, it is convenient to consider how many and which 
quantities may be derived from ionospheric measurements. Records of iono- 
spheric conditions are almost universally given in the form of graphs of virtual 
height plotted against frequency; these are obtained either by automatic 
equipment such as that developed by Berkner, Wells and Seaton (1936) at the 
Carnegie Institute, or by a more modest form of manually operated equipment, 
where financial resources are not readily available or the proposed experimental 
work is not of a very permament nature. From such records may be obtained: 


(i) the maximum electron density N in the layer, which is directly proportional 
to the square of the critical frequency for the ordinary ray ; 


(ii) the true height, hy, of maximum density, and semi-thickness 7 of the 
layer, which may be obtained without much difficulty either by the 
method developed by Booker and Seaton (1940) or by a more recent 
method developed by Pierce (1947). 

At this stage it is interesting to remark how investigators have, up to the 
present, insisted on jumping to conclusions, basing their reasoning entirely 
on critical frequencies and virtual heights. In dealing with a problem as complex 
as the ionosphere, it is an obvious advantage to rely as far as possible on 
experimental results and thus to eliminate, in part at least, some of the little 
known parameters. 

The photographic records used for purposes of testing this theory and 
gaining some information of conditions in the F, region were taken in this 
laboratory with manually operated equipment assembled by the writers during 
the war. A pulse transmitter of continuously variable frequency covered a range 
from 1:8 to 13-9 Mc/s.; its mean power output was about 200 watts during pulses. 
Records were made on 35 mm. film and scaled through a microfilm reader giving 
an approximate magnification of 10. By the help of frequency marks made on 


Ionized Layer in Non-Isothermal Atmosphere 435 
the record at regular intervals of 0-2 or 0-5 Mc/s. by means of a height calibrating 
oscillator working at 1,500 c/s., frequencies could be read to 0-01 Mc/s. and heights 
could be reproduced to 2 km. 

A series of photographic records was taken, one on each hour on alternate. 
days between Ist November 1945 and 28th February 1946, the summer season 
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in South Africa; this choice of intervals was dictated by the shortage of observing 
staff available. All the records were completely scaled according to the method 
of Booker and Seaton (1940); heights were read at the 8 frequencies suggested 
by Booker and Seaton, and coupled in four pairs, each of them giving the true 
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height of maximum density, 4, and semi-thickness, 7, of the layer. The mean 
of the four values was calculated for each quantity. Where necessary, corrections. 
due to retardation by lower layers were applied (Booker and Seaton 1940). 
. Since a significant correlation was found between magnetically disturbed and 
ionospherically disturbed days, all days with magnetic character figure 2 were 
eliminated in taking the monthly mean values of N, hy and 7. As some large 
deviations still remained, the more outstanding ones were removed, using the 
well established test that deviations larger than three times the probable deviation 
are unlikely to be due to chance alone. With the adjusted monthly mean values 
for N, hy and 7 thus obtained, it was possible to draw the Booker and Seaton 
parabolic distribution of (electron density, true height) for each hour of the 
average day for a given month. 

It has been found convenient to represent all the information about the 
behaviour of the ionosphere on the average day for each month by a graph showing 
lines of constant ” on a coordinate plane of time and height. This plot is called 
an ‘electron density map’: the map for the February data is shown in Figure 6 
as an example. It serves to illustrate the important point that dn/dt is in tact 
nearly zero during the daytime in parts of the F, layer below the electron density 
maximum. 


§3. APPLICATION OF THE THEORY TO EXPERIMENTAL RESULTS 
3 (i). The Method of applying the Theory 


‘The application to experimental results of a theory with equations as complex 
as those derived in § 1 of this paper is a problem which must present some difficulty. 
From the theory there emerge the three useful equations (18), (22) and (24). 
Considering C as known, we see that these equations contain no less than four 
unknowns: B, y, Ty) and F. B may be eliminated by dividing (18) by (22), 
thus obtaining (28), but this leaves only two equations with three unknowns,, 
Vy bg atid E, 

At first sight it might seem that by taking two values of and the corresponding 
two values of n from the experimental parabola, we could obtain two numerical 
equivalents of (28) which, together with (24), would allow solution for all three 
unknowns. ‘This is a fallacy, however, because by dividing equations (18) and 
(22), the absolute value of N has been eliminated and the parabola now yields. 
only two quantities. 

It is therefore necessary to assume a numerical value for one of the three 
remaining unknowns, the one chosen being that which can best be assessed. 
The temperature 7) is obviously the one to be chosen, as it is not likely to be 
outside the limits of 350°k. and 700°k., i.e. it varies through no more than a 
factor of two. We shall not, however, assume a fixed value of Ty, but apply 
the theory for a number of values of T, between these two limits, deferring 
discussion of its actual value until later, 

Algebraic solution of equations (24) and (28) is impracticable as a routine 
method, and resort is therefore made to a graphical method, which can be applied 
quickly to most conditions met with in the F, region. It should be remembered 
that the final equations are only true if dn/dt =, a condition which is well fulfilled 
for the F, region below its maximum from 08h. to 16h. South African Standard 
Time, as previously pointed out. 
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3 (ii). Graphical Solution for F and y 


The value of n for substitution in equation (28) may be read at any height 
below hy, where the electron density is greater than 0-3N. ‘These limits are those 
beyond which the curve differs appreciably from parabolic form, and between 
which all the virtual heights are originally scaled. 

In order to construct a family of curves of (n/N, y) for various values of F' 
from equation (28), the value of (k—h»)/T) must be settled. A value of h 
agreeing with the conditions outlined in the previous paragraph is 280 km., as 


. n : . hh i j 
Figure 7. Curves of N against y for ag =e from equation (30). 


may be verified by an examination of the electron density map. Assuming that 
T, will be of the order of 400°K., we find 

h—h, 280-200 1 

Sie ore oe Se (29) 
When we recall that C=16:2°c. km“, and substitute in (28), dropping the 
subscript b, this simplifies to 


n\2 a y — (1+ 16+2/y) ( y = (1+ /16-2) 
() =(1+3) F(1+7g5) 
; y 1 y —16-2/y 
x exp (1 = ts) (| a5 4 ot pg ay ee (30) 


Choosing values of y which give simple numerical values to the indices 
i= 1-62, 3-24,.4°36, 6:49, 8-10, 9-7), and using equations (15) and (21) for y=0, 
various values of F may be substituted to give n/N greater than 0-3. In this way 
the family of curves shown in Figure 7 was drawn. 
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Equation (24) may be written, with the usual value of C, 


“iene z auc +75) fee -1]. — (31) 


With the values of y used before a similar family of curves may be drawn, and is 
shown in Figure 8. 


These curves exhibit a number of interesting points: 


(1) The occurrence of maxima in the curves of n/N for F=5 x10 and 
F=10" for values of y equal to 1-2 and 8 deg.km™ respectively 
indicate that under these conditions m becomes equal to N, i.e. the 
maximum electron density passes through the height 4 for which 
(h—200)/T, =3. 

As was illustrated by Figure 5, an increase in F' corresponds to a 
lowering of the height, and the curves for values of F greater than 10-4 
refer entirely to values of n above hy, and therefore cannot be used 
for practical purposes. 


(ii) In Figure 8 the curves below (Ay, —200)/T)=0-2 again refer to conditions 
when the maximum is below the selected height h, and therefore have 
no practical value. 
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The simultaneous graphical solution of equations (30) and (31) may now be 
accomplished, if values taken both from Figure 7 and from Figure 8 are plotted 
on the same diagram, using as coordinates the variables F and y or, preferably, 
log F andy. ‘This has been done in Figure 9, where lines of constant n/N from 
Figure 7 are drawn in full and those of constant (hy —200)/T, from Figure 8 are 
shown broken. 
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The process of evaluating y and F for a given experimental parabola now 
resolves itself into the following series of simple operations: 


(i) select a reasonable value of Ty; 
(ii) estimate the value of h by substituting this value of 7) in the equation 
(h—200)/Ty=3; 
(iii) from the parabola, find 7 at this value of A, and also N; 
(iv) hence calculate n/N; 
(v) from the parabola read off hy, ; 
(vi) evaluate (h,,—200)/T, using the value of Ty chosen in (i) ; 
(vii) on Figure 9 locate the point which satisfies simultaneously (iv) and (vi), 
and read off the values of y and log F. 


Log £ 


0-6 


Lines of constant 47 
Ruy~No 


7 (°c.km-') 


Figure 9. Graph for the evaluation of y and F from experimenta! data. 


The whole procedure is repeated for each value of 7 selected. 

The curve n/N =1-0 represents the maximum electron density, and thus the 
part of the graph above this is not shown, as it has no use. It coincides with the 
line (Ay, —200)/T, =0-2, which also represents the electron density maximum. 

The angle between the two families of curves increases as we estimate 1 
further below the maximum, and therefore values of y and F estimated from 
values of scaled near the maximum are not very accurate. ‘This could be 
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avoided by choosing the fraction (h—200)/7T)A+, but this would introduce other 
difficulties by putting 7 below 0-3. in the normal applications. 

The two families of curves tend to run parallel for values of y greater than 
15 deg.km~?, and this introduces uncertainty in values estimated in this region. 
Since for y=15 deg.km~! and T7)=250" k. (a very low estimate), the temperature 
at 350 km. would be 2,500°k., which is more than anyone has hitherto dared to 
‘suggest, this is no disadvantage. Godfrey and Price (1937) have shown by 
considerations of radiation equilibrium that this temperature cannot rise above 
3;300°R, 


3 (iii). Treatment of Experimental Results 


The process just developed was applied to the experimental results for 
November 1945—February 1946, described in §2 of this paper. The value of T, 
‘being unknown, the estimation of y and F was repeated for each parabola for a 
series of values of 7), varying in steps of 50°x. from 250 to 600°x. In many 
cases the higher values of 7) gave values of y well off the graph, and these, in 
accordance with the previous discussion, have not been taken into account. 

To illustrate the results of this treatment and the following discussion, Table 2 
as given. It shows, in the first two rows of each section, the figures for y and F 
for November. Similar tables were prepared for the other three months. 

The values of y obtained are practically all above 5 deg.km~ which points 
straight away to temperatures of the order of 1,000°K. in the F, layer. Further 
discussion on this point is, however, premature, because we have no means of 
selecting any particular value of T, as yet. 

Some indication of the correct value of Ty and y may be obtained by 
performing two tests. 


Lest £. 
From equation (17) we find 


_. Ccosx 
TTR ee: On eg ee ( 


As the maintenance of the layer is presumably due to ionization of a particular 
species of atom or molecule, 8 should be constant. Cosy may be calculated 
from the Nautical Almanac and we have a value of F for each Ty. However, 
ny will also vary with Ty. During the preliminary work this was calculated from 
a probable vaiue at 100 km. on the assumptions that there the temperature 
remained unchanged at 360°K. throughout the day and that there was a linear 
temperature gradient between the 100 and 200 km. levels. Later a more 
detailed examination of the problem gave what is believed to be a better value 
for the 200 km. level, mj =1-:5 x 10!" cm-3 when 7) =360°x., but did not indicate 
the variation of this quantity with 7). For the present purposes this quantity 
has been taken as a datum and modified for other values of Ty by assuming it 
to vary in the same proportion as the earlier estimate. The values so obtained 
are shown in ‘Table 3. 

In this way the values of 8 corresponding to each T, at each hour were 
evaluated, and those for November are given in the third row of each section 
of Vable 2: 
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Table 2. Experimental values of y and F and derived values of 8 and B/nj. 
November 1945 


S.AcS: By (ea. | B5250 300 350 400 450 500 550 600 
Time 

y 11-7 10:7 9-2 8:5 8-0 9°5 11 — 

8h. 100F 3:94 4:34 4:68 SoD Soo) 7:50 -9-33 — 

cos y=0°529 Bx10 10-6 6:04 3:63 2:41 1-82 1:17 O:78 — 

Bing 4940 2740 1440 901 675 529 396 — 


y 11-4 9-9 8-9 PS UES 6°5 10 — 
Qh: 100F SV 4702) 8440-5 © 458725270)" 95°S90 87 

cos ¥=0-701 Bx 10e >) -15°0 S052 600 So 4a 0.248 . 1-97e el10 5 
Bin, COLORS 120 OOOReLOSOME OLS 565 461 — 


y 11325 ibe ah ABIES: ORR - Oph hihes 13-0 = 

10h. 100F A300 501 Say See @dds Wey Oeil 

cos x¥=0°836 px 10%=— 13-1 B27 eee Oo a8. 2c O0 ae e768 (loDy = 
B/no St 70m 42405 2370215505 1160942 dS — 


y ass) OWS) 883 ISS Fel 7-0 55) == 
11h. 100F 075 3°30) e301 AGL 4b He 7il SNS) = 

<os x=0-927 pol 0e S228 137 8:25 oN hoi A ICTS 
B/no (NCO) GRD SSS AD ANSTO AAD 890 — 


y 10:2 O53 8-0 ie OLO 6:7 Soe) Ws 
12h. 100F 2-69 20S :145 = 3°36, 3°70" 4-20 §4:98) 925-07 8:32 

cos x=0-966 Bx 101 2228:3 TESS) 9220 ee 618.284 O50 = eso h BX ANS 
Bing 47200 8470" 4520- 2860 21107. 1510 1110 ~ 920 


y 8°5 8:2 Wl 6:0 6:0 Byes 4-0 9 
ish: 100F Datel Zoi 22 19 03:82 422 eel ee 4, 

cos ¥=0-952 B><e10% 35-8 “4-17-6. 10-4 e/a Oi Beet Baily and) 
B/no 15600 8070 4210 2600 1920 1390 993 812 


y 9°8 8-9 7°5 ee} (HPI! 5*4 SOMOS O 
14h. 100F 20510 0a 29 e873 13. 64440 4-09 3-03 

cos ¥=0°884 pile 25-6 14-4 SI 270) aoegiil BE) So hah) 
B/ng 15300 =. 69908) 36802520. 91690 41210 895 720 


y 7°8 Ts, 6°3 DOMMES 3°8 8 — 
15h. 100F 2e 257 9 wes 09 355s 80 Be 4 Oa 87595 

cos ¥=0:769 ped tae 270 13-6 7-9o=— 5-19 4:01 ls) (hese) ee 
B/ng 11200 5770 3000 1900 1370 960 726 — 


y 6°5 6:1 5c) 4°8 3:6 5 12 — 
16h. 100F 2:01 Zoi. 2°98 SL ono 5°62 10 — 
cos y=0°613 p10)" 24-0 1 G5 Oct /ar oe SS) 1-81 0:34) —— 


B/no 9020 4520— 2310 1440 1040 718 689 — 


Table 3. Values of nj for various 75 
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Test II. 

The second test depends on the evaluation of the constant B, previously 
eliminated by division of equations (18) and (22). B is defined by equation (16) 
and is also given by a rearrangement of equation (22). From its definition it is. 
seen that B is not a constant, since it contains /, which is a function of Ty. The 
quantity 


Bom ile ' 
ee from equation (16) 
N? exp (1+ y/C) 


ws he frog Cc anions 2c) ame eee (33) 


= (i+y7/0) 
nf F(1 in zr) é 


should however be constant at all times and for all months. The values for 
November are shown in the fourth row of each section of Table 2. 


3 (iv). Discussion of the Tests 


The correct value of B is very doubtful. Table 2 and the corresponding 
tables for the other three months have been drawn up to cover temperatures T) 
from 250° x. to 600°K., which range includes most of the rough estimates found _ 
in the literature, and, if this range is accepted as the most probable for To, the 
correct value of 8 should appear in the third row of each section. An examination 
of the tables shows that values of 8 in the neighbourhood of 3 x 10-1” cm? fulfil 
this condition. Table 4 shows the values of 7), interpolated to the nearest 10° k. 
from the tables, which correspond to 8 =3-2 x 10-1!” cm?. 


Table 4. Values of T, (°K.) corresponding to B =3-2 x 107" 


SINS 

Time November December January February 

(hours) 
8 370 320 320 320 
9 410 380 350 360 
10 410 410 470 390 
dal 480 540 410 430 
2 500 490 470 520 
13 540 460 430 450 
14 500 450 420 440 
15 490 430 390 410 
16 450 390 330 390 


In theory, if a value of B/nj is chosen to give the same 7 as the value 
B=3-2 x 10-” cm? at a certain time, this same value of B/ng should correspond 
to the identical set of values of 7, shown in Table 4. However, B contains « 
and J,,, which may vary with season, sunspot activity, etc., respectively. Hence 
we shall not expect exact: agreement in practice. ‘To find the value of B/ng 
which best corresponds to 8 =3:2 x 10-!" cm?, a figure was estimated to correspond 
to each of the noon values of 7, in Table 4, and the mean of the four was taken. 
This gave the best value of B/nj as 910 cm™*; the agreement was not good, 
the largest deviation from the mean being about 50%. Considering the 
relatively large percentage change in B/nj corresponding to a change of only 
50°. in Ty, however, this may be regarded as satisfactory. The J) values 
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corresponding to B/nj=910 cm-’, interpolated for the last row of Table 2 and 
the corresponding ones for the other three months, are shown in Table 5. 


Table 5. Values of Ty (°K.) corresponding to B/nj,=910 cm? 


S.A.S. 

Time November December January February 

(hours) 
8 400 330 300 320 
9 430 330 330 400 
10 510 350 350, 450 
ill 550 440 350 450 
12 600 410 370 550 
13 570 410 350 560 
14 550 420 DN) 550 
15 510 390 340 550 
16 470 380 310 460 


The agreement between the two tests is very pleasing, as they rest on 
independent grounds; therefore it does not seem out of place to take the mean 
values from these two tables the best representative values for further discussion. 
These are given in Table 6. 


Table 6. Mean values of Ty (°K.) 


S.A.S. 

Time November December January February 

(hours) 
8 385 325 310 320 
9 420 B55 340 380 
10 460 380 410 420 
11 515 490 380 465 
12 550 450 420 535 
13 555 435 390 505 
14 525 435 385 495 
05) 500 410 365 480 
16 460 380 320 425 


The absolute values of these T,’s depend, of course, on the arbitrary choice 
of Bas3-2x10-"cm?. The effect of changing this was investigated by estimating 
new values of JT, by the two tests with B=1-6 x 10-1” cm? and 8-0 x 10-2? cm? 
respectively, and the corresponding figures for B/nj. This resulted in an 
increase of the estimated 7, by about 100°x. in the first case and a decrease by 
the same amount in the second case, but it did not affect the form of the variation 
of 7, with time of day, as may be easily verified by an examination of Table 2. 

Previous estimates have been made of the mean temperature at 200 km., 
but they vary widely in order of magnitude. Martyn and Pulley (1936), from 
considerations of electron collision frequencies and of auroral spectra, have 
tentatively adopted the value of 1,000°K. at 200 km. in the daytime; Hulburt 
(1939b) has assumed the value 360° k. above 100 km. for lack of better evidence; 
while Vassy and Vassy (1942), assuming complete dissociation of both oxygen 
and nitrogen, obtained a temperature gradient which gives T) at 200 km. as 
600°K. in summer. ‘The values in Table 6, having a mean of about 430° K., are 
well in accord with these estimates. 

It may be seen that Ty rises towards the middle of the day and falls off with 
some delay in the afternoon. ‘This is true no matter what f is chosen to be. 

2 H-2 
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It shows that the heating of the atmosphere does take place even at a height of 
200 km., the temperature rising by about 30° of its value at 8 hours. 

Estimates have been made of the ‘temperature of the F, region’ during 
summer daytime: 


(i) Appleton and Naismith (1935), from considerations of the variation in N 
between summer and winter noon, estimate that the temperature 
cannot be less than 1,200° k. in summer. 


(ii) Fuchs (1936) has concluded from estimated scale heights that the 
temperature of the F, region may rise as high as 1,900°K. 

(iii) Martyn and Pulley (1936) confirm Appleton’s result of a temperature of 
1,200° k. 

(iv) Vassy and Vassy (1942) assuming complete dissociation of O, and Ng, 
obtain a summer temperature of 1,080° k. from the estimated scale 
height of 70 km. given by Appleton. 

(v) Godfrey and Price (1937), by a consideration of the equilibrium between 
absorption and secondary radiation, conclude that the temperature 
is of the order of 1,000°k., and cannot exceed 3,300°K. 

For comparison, the temperature at 300 km. was calculated from the T)’s of 
‘Table 6 and the corresponding y’s. ‘The individual variations of these figures 
for Tso) are probably not significant of actual conditions, and therefore only the 
means are included in Table 7. 


Table 7. Mean Temperatures at 200 and 300 km. 


Month November December January February 
Mean Ty (° kK.) 485 408 369 447 
Mean T%o9 (° K.) 1097 1376 1392 1305 


It will be observed that the mean values for TJ, are smaller for the midsummer 
months December and January than for November and February, while the 
reverse holds for temperatures at 300 km. ‘This gives valuable support to the 
expansion theory, for, near the solstice, when heating is greatest, the atmosphere 
at 200 km. becomes heated in the middle part of the day, and is forced upwards 
by expansion, leaving cooler gas to take its place. ‘The heated mass of gas would 
produce greater absorption of the solar radiation, thus preventing it from reaching 
the 200 km. level, where it would normally cause more heating. 

As a last interesting point, we may make an estimate of the intensity of active 
incident radiation [,, above the earth’s atmosphere. From equation (16) 


(ee (=) “3 Se (34) 


We have assigned to B/nj the value 910 cm from the results of test I. The 
value of « is open to discussion, and seems to vary with ion density and height, 
according to Wells and Shapley (1946); a reasonable mean seems to be 
5 x 10-1 cm sec-1; B has been fixed at 3-2 x 10-1”? cm2; and the value of the 
ionization potential w is given by Bacher and Goudsmit (1932) as 13-55 ev. for 
atomic oxygen and 14-48 ev. for atomic nitrogen. As a mean 14 ey. has been 
used. ‘Thus we find 


J. =0°313 ergvom=atiséca) 
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- Hulburt (1939b) has calculated that the total black body radiation of energy per 
quantum greater than 13-7 ey. received from the sun is 8 x 10-3 erg.cm~?sec“!, 
which is much lower than the above estimate. 

Substituting back the value 8 x 10-% in equation (34), we find B/nj =23-2, 
and from ‘Table 2 this will be seen to correspond to 7, at least as high as 2,000° k. 
This is ridiculous, and we can only conclude that the sun does not behave as a 


black body in this spectral region. A similar conclusion has been drawn by 
other workers (Bradbury 1938, Saha 1937, McNish 1937, Hunter 1943). 


VAG ONCE WU STOIN'S 

The introduction of a linear temperature gradient into the theory of layer 
formation yields equations which may be applied with success for data for the 
F, layer, presented in the form of parabolae representing the distribution of ion 
density with height. ‘The values of the temperatures at 200 and 300 km. obtained 
are in accordance with previous estimates, and are of a more quantitative nature. 

The main weakness seems to lie in the assumption regarding the absorption 
coefficient 8. In the first place its value is not fixed by the theory, and in the 
second place it varies with the frequency of the radiation absorbed, and an 
assumption of a mean value is not necessarily justified. A general treatment of 
this problem has been given by Chapman (1939). Again, the manner in which 
the recombination coefficient « varies with the temperature and pressure is not 
wellestablished. Finally the theory is limited in its application by the assumption 
that dn/dt=0 and by its neglect of the earth’s curvature. Re-examination of 
these points seems to be the next step in its development. 

It should be pointed out that the assumption of complete dissociation of 
oxygen and nitrogen above 100 km. is not a weakness, since a numerical 
recalculation is all that is required to apply the theory to an atmosphere containing 
constituents of known molecular weight. 
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A New Method of Measurement of the Variation with 
Wavelength of the Refractive Index and Absolute 
Stress Optical Coefficients of Amorphous Solids 


By J. E. H. BRAYBON 
University College, Southampton 


Communicated by S. Weintroub; MS. received 9th Fanuary 1950 


ABSTRACT. White light is passed through a Michelson interferometer in one of whose 
arms a thin lamina of the substance in the shape of a tension member is supported by a 
loading frame. The emergent beam, examined by a spectrometer, shows a continuous 
‘spectrum crossed by interference fringes. For the unstressed lamina the position of the 
fringes gives the variation of refractive index with wavelength. For the stressed lamina 
the behaviour of the fringes gives the stress optical coefficients. No modification of the 
interferometer, and only three sets of observations on the fringes, are required. For the 
calculation of absolute values of the refractive index the values of two wavelengths must 
be known. 

Full details of the method are given, and photographs of tests on Columbia Resin 39 are 
included. These show that CR39 becomes positively birefringent under tension. 


§1. INTRODUCTION 

HE method, which uses an unmodified Michelson interferometer, has been 

specially developed in order that measurements over a wide range of wave- 

length of both the refractive index and stress optical coefficients can be made 
on the same specimen. ‘This is of importance in the case of a plastic, where the 
optical properties may vary with age and from sample to sample, and thus, unless 
the same specimen is used throughout, the refractive index cannot be correlated 
with the stress birefringence. The material used in checking the method was a 
condensed allyl ester, Columbia Resin 39 (CR39). In one arm of the interfero- 
meter, a thin lamina of the material (1 to 2mm. thick), in the form of a tension 
member of uniform cross section, is held by a frame which may be used to apply 
a homogeneous tensile stress to the specimen. White light is passed through 
the system and the emergent beam is examined by a spectrometer. ‘The con- 
tinuous spectrum is observed to be crossed by a system of interference bands. 
From this system for the unstressed lamina, the variation of refractive index with 
wavelength is obtained; from the behaviour of the system when the lamina is 
stressed the stress optical coefficients are deduced. The theory of the method 
is given in detail. 


§2. DETERMINATION OF REFRACTIVE INDEX DISPERSION CURVE 
(1) Experimental Details 


The material, in the form of a plane parallel thin lamina, is inserted in one 
arm of a Michelson interferometer. Light from a ‘ Pointolite’ source P (Figure 1) 
is focused on an adjustable circular diaphragm I, collimated by the lens L and 
passes into the interferometer M. Emerging from M, it is focused as a vertical 
line by the cylindrical lens C on to the slit of a Hilger constant deviation spectro- 
meter H, or a spectrograph. 
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If great care is taken to align all the components strictly along the optical axis 
of the system, then, for a suitable setting of the tilt of the interferometer mirror 
Mg, the resultant spectrum is observed to be crossed by vertical dark interference 
fringes. The relation between fringe wavelength A, refractive index », and 
thickness ¢, of the lamina may be expressed as 


Be DA her LeMay ee hl eee (1) 


where D is the optical path difference between the two interfering beams in the 
absence of the specimen and 2, the fringe order, is an integer. 


m, 


See ee 


Figure 1. Optical system. 


Due to the dispersion of » witha the fringes are not equidistant in wave number 
but are approximately symmetrical about a wavelength which is dependent upon 
D, i.e. upon the setting of the movable mirror, m,, of the interferometer. ‘The 
photograph, Figure 2 (Plate I), obtained with a lamina of CR39, thickness 
1-68 mm., shows the variation in position of this central fringe with alteration of 
D by lateral movement of m,. Further, replacing the lens C by a spherical 
lens, the diaphragm I by a vertical slit, and the source by a line filament, circular 
fringes are obtained and, in the region of the central fringe, it may be confirmed 
that corresponding fringes on either side of the central fringe are of the same order. 
Also, from equation (1), it is apparent that the fringe order rises to a maximum 
at the central fringe. 


(ii) Calculation of the (1, A) Dispersion Curve 

The absolute value of the fringe order cannot be determined directly without 
knowing the other variables in equation (1), but by plotting the fringe order 
(arbitrary) against A, a curve of the form shown in Figure 3 can be obtained. From 
this curve the dispersion curve of j against \ may be obtained as follows: Let the 
wavelengths of two fringes, on opposite sides of the maximum, but having the 
same fringe order, x (absolute value unknown), be Ay and A,’._ If the wavelength 
of any fringe of order (« +m) isA,,, and jp, 4’ and yu, are the values of the refractive 
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index at Ay, Ap’ and A,, respectively, then, from equation (1), by elimination of x 
and D: 

Ho — Man = (to — Ho’ )(Am — Ao) /(Ao’ = Ao) + Am/2E- sev (2) 
If yy and pip’ are known, j,, is given by the above equation. However, it is not 
necessary in practice to calculate .,, directly from equation (2), as the (1, A) 
dispersion curve can be obtained by a simple graphical construction. 

In the absence of the mA,,/2t term, equation (2) represents the broken st.aight 
line of Figure 3 (5), through the points (Ag, 49); (Ao » Ho )» Which is easily constructed. 
By adding to the ordinate of any point on this line the corresponding value of 
mA,,/2t the corresponding point on the dispersion curve is at once obtained. 

It should be noted that this interpretation of equation (2) confirms that the 
fringe order rises to a maximum, since, in order to obtain the normal dispersion 
curve, the term mA,,/2t must be positive in the range Ay to Ay’, and negative outside 
this range. 

For this method, either 4) and jy’ must be determined by other means, or D 
and one value of refractive index must be known. With the interferometer used 


Fringe Order nr 


Wavelength A 


Refractive Index “4 


Wavelength A 


Figure 3. (a) Variation of fringe order (7) with wavelength A. 
(6) Variation of refractive index (1) with wavelength A. 


it was more convenient to adopt the former course. Furthermore, unless the 
lamina itself is absolutely plane parallel, the interferometer, when the lamina is. 
removed, will not give fringes without an alteration in the tilt setting of the mirrors 
in addition to a change in D. 'This completely precludes a determination of 
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D since an adjustment of the mirror tilt in turn affects D. Because of this yy and. 
fy Were determined on two other specimens, both in the shape of 60° prisms. 
Experience had shown that in the case of CR39 the use of other samples is. 
permissible as its optical properties are quite constant. Such a course, however, 
would not be suitable for a plastic such as Catalin, which is very susceptible to time 
effects, and the determination of «4 must be made in this case directly by a refracto- 
meter or by a method involving the measurement of the Brewsterian angle 
(Pfund 1941). Thefinal dispersion curve obtained agrees very closely with readings. 
made using the 60° prisms at various other points over the visible range. ‘The 
lamina used was supplied by the makers from their normal range, and, apart. 
from cutting it into the form of a link for the purpose of stressing it as described. 
below, it required no treatment. 


$3. MEASUREMENT OF THE STRESS OPTICAL COEFFICIENTS 
(i) Experimental Details 


The lamina used must be in the form of a tension member which is subjected. 
to a homogeneous tensile stress by means of a simple bell crank loading device. 
On loading as indicated in Figure 1 the fringe system described above is displaced 
outwards from the central fringe and the fringes over a wide range of wavelength 
disappear as the load increases. However, fringes are again visible if a polarizing 
prism is placed in the beam of light entering M, and the prism turned so that the 
longitudinal axis of the specimen (and therefore the direction of the stress) is. 
parallel or perpendicular to the direction of the electric vector. 

These effects occur because the specimen becomes optically anisotropic under 
load with its axis of symmetry along the line of stress. 

The birefringence induced by the stress results in the fringes observed in light 
polarized with its electric vector at right angles to the tension being displaced 
by a different amount from those observed in light polarized with its electric 
vector parallel to the tension. Fora suitable stress, over a comparatively large 
wavelength range the minima of the fringes obtained in the one position of the 
polarizer correspond to the maxima obtained in the other position. Thus, in 
the absence of the polarizer, no fringes are visible. If the load is further increased 
the fringes again become visible, indicating that one system has continued to move 
relative to the other. 

This is illustrated in the photograph, Figure 4 (Plate II), which was obtained 
by placing a Rochon double image polarizing prism in the collimated beam entering: 
the interferometer, and a spherical lens in place of C for focusing the resultant 
two beams from M, polarized at right angles to each other, upon the spectrometer 
slit as two separate point images. A double fringe system is then produced. 
In the absence of any stress the two sets of fringes are in alignment. ‘The fringes. 
move away from the central maximum and move relatively to each other as the 
load is applied. 

It is of interest to note that the spectrogram taken one minute after the removal 
of the load shows that the fringe systems have not yet returned completely to 
coincidence. This indicates some residual birefringence due to creep (Braybon 
1949). That there was also a very small amount of birefringence present before 
the load was applied may be seen from close examination of the spectrogram taken 
prior to loading. This was due to the immediate past history of the specimen, 
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(ii) Calculation of Stress Optical Coefficients 


Photoelastic properties of an amorphous solid are characterized by two stress 
optical coefficients B, and B,, which at any given wavelength may be defined — 
(Mueller 1935) by 


(Hy — pH) =Sp,=B,P,, wa ae (3) 
(4,—) =du, =BP,, 


where P, is the pressure in the z direction. | 


u is the refractive index of the unstressed medium and yp, and j, are the re- 


sultant indices of refraction in the stressed substance for light with its electric | 


vector perpendicular and parallel to the direction of pressure. The sign of P, 


is thus negative for a tension. By subtraction: 


B=B,—B,=(u,—p,)/P, = hues (4) 


where B is the relative stress optical coefficient. Most determinations are of the 


stress birefringence and merely give the constant B. 
At any given wavelength, the values of 5u,, and du, may be deduced, for 
‘differentiation of equation (1) keeping A constant gives 


See (A128) Oh th. a (5) 


The changes in fringe order, 5n,,5n,, for light polarized with its electric vector 
respectively perpendicular and parallel to the stress P., are found by comparing 
(fringe-order, A) curves obtained as above with those obtained in the unstressed 
‘condition. ‘The corresponding refractive index changes are therefore: 


Se = —(A/28)8n, = BP. | 
dp, = — (A/2t)in, = BSP., 
and the resultant birefringence Aw is given by 
At= (Ou, —du,)='— (6n,=0n,)A/2t=—- BPA Bee (7) 


It should be noted that for a tension P, is negative; for CR39 it was found that 
én, and dn, are both positive and that 5”, >5n,. Thus positive values of B,, B, 
and B are obtained, indicating that, as with the majority of other substances 
‘showing stress birefringent properties, CR39 becomes positively birefringent 
under tension. 

Furthermore, if the total load applied to the specimen is L, equations (6) and 
(7) may be rewritten, with 1/2¢P, replaced by b/2L, where b is the breadth of the 
specimen, a quantity which may be determined with less error than t, which is 
of the order of only a millimetre. 

The application of a load will result in a change in the dimensions of the 
specimen, i.e. in b and ¢t. For the low values of load employed, however, it was 
considered unnecessary to correct for this variation which for the load used is 
of the order of 0-05. 


§4. DISCUSSION 
The method for the determination of stress optical coefficients bears some 
resemblance to that of Mach (1872) and of Twyman and Perry (1922). However, 
in the first of these two methods a Jamin interferometer was used and the method 
depended upon the use of a compensator block of the same material and upon a 


| 
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second measurement with specimen and compensator immersed in a fluid of 
known refractive index. In the second method a Michelson interferometer was 
used but the specimen was subjected to flexure and examined with monochromatic 
light. 

The advantages of the present method are that only three sets of readings are 
required; measurements are made for small values of stress, reducing creep 
effects to a minimum; the specimen need not be very thin and fragile and its 
dimensions can be determined accurately. Estimation of the fringe shift for 
values of <A, is somewhat difficult owing to the rapid rate of decrease of fringe 
order with decreasing wavelength, and results for the range may be somewhat less 
accurate. Obviously the method requires that the straight line through the points 
(Xp, feo) and (Ap’, 49’) be accurately determined and carefully drawn. 

The method has been used successfully on CR39 over a wide range of wave- 
length. The work is being extended to other plastics and full details of the 
results will be published later. 


ACKNOWLEDGMENTS 
The author wishes to thank Professor A. M. Taylor for his advice and 
encouragement, and Mr. S. Weintroub for his assistance in the preparation 
of this paper. 
REFERENCES 


BrayBon, J. E. H., 1949, Nature Lond., 163, 689. 

Macu, E., 1872, Ann. Phys., Lpz., Ser. II., 146, 316. 

MUELLER, H., 1935, Physics, 6, 179. 

Prunp, A. H., 1941, ¥. Opt. Soc. Amer., 31, 679. 

TwyMan, F. and Perry, J. W., 1922, Proc. Phys. Soc. Lond., 34, 151. 


Vibrations of Free Elliptical Plates 
By MARY D. WALLER* 


MS. received 16th January 1950 


ABSTRACT. The systematic observations which have been made on the normal 
vibrations of free elliptical plates confirm Chladni’s results. "The nodal systems may be 
divided into the same four classes of symmetry as those of rectangular plates. 

The correspondence which exists between the vibrations of elliptical, rectangular and 
circular plates is traced. Combined modes of vibration do not occur on elliptical plates as 
they occasionally do on rectangular plates. 


Si, INTRODUCTION 

HE only extensive experimental study of vibrating elliptical plates would 

appear to be that of Chladni. The 23 drawings given in Die Akustik, 

1802, concern a plate whose axes are in the ratio of about 3/2 and the nodal 
systems are arranged systematically, as in the photographic records of the present 
paper. The text, §§ 142-157 includes data regarding the vibration frequencies 
of many other plates which may be converted from the notation of the chromatic 
scale by using the Zamminer table given in Rayleigh’s Theory of Sound, I, 
1894. 
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Pavlik (1937) has used the Rayleigh—Ritz approximate method (Rayleigh 1894, 
1911, Ritz 1909) to calculate the first seven natural frequencies of a rather wide | 
stainless steel plate (ratio of axes 1:29) and to compare these with observation. — 
The agreement is good for the lower frequencies, but for the sixth and last pair 
of frequencies the values given by theory exceed the observed values by 44° and — 
52% respectively. ‘The theory, which involves laborious calculations, has not 
been tested on narrower ellipses. | 

It is evident that our information regarding the vibrations of elliptical plates | 
must depend principally upon experiment, and it is desirable that photographs. 
of the nodal designs produced should be placed on record. 

The present work forms part of a more general study of the vibrations of plates _ 
which has not as yet reached completion. Details regarding the experimental 
technique have been given in previous publications (Waller 1949, and references _ 
therein). | 

The subject is approached through considerations of symmetry which make _ 
it easy to arrive at a general understanding of vibrations of plates of any given | 
geometrical shape. 


§2. SYMMETRY AND CLASSIFICATION OF THE NORMAL 
NODAL SYSTEMS 

The elliptical plate possesses two mirror lines of symmetry, the major and the 
minor axis, and the rotational symmetry is 180°. Lines of symmetry divide the 
nodal pattern so that one half is the mirror image of the other. The fact that the 
lines of symmetry may be nodal or antinodal may be appreciated by examining the 
first nodal system of the first two rows of Figure 1 or 2 (see Plate). In the 1]1 
system of the second row the lines are nodal and the displacements of corre- 
sponding points of the surface on either side of them are equal and out of phase. 
In the 0|2 system of the first row the lines are antinodal, that is to say they pass 
through the points of greatest displacement, and the displacements of corre- 
sponding points on either side of them are equal but in phase. 

It is evident that all possible nodal systems can be divided into four classes 
according as (i) the centre is antinodal, (ii) the minor axis is nodal, (111) the major 
axis is nodal, (iv) both the axes are nodal. 

‘These classes correspond with those of rectangular plates where the two lines 
of symmetry are the two medians (Waller 1949). 

All the nodal designs can be brought into self-coincidence by a rotation of 


180° about the axis passing through the centre of, and perpendicular to the plane 
of, the plate. 


§3. THE NORMAL NODAL SYSTEMS 

Photographic records of the nodal designs obtained on a narrow plate (ratio 
of axes 2|1) and a wide plate (ratio of axes 5|4) are shown in Figures 1 and 2 
respectively. 

‘The arrangement of the nodal systems is the same as Chladni’s and it corre- 
sponds also with that of rectangular plates (Waller 1949). 

Fach column, as indicated by the number above it, contains systems with a 
constant number of lines running in the direction of the shorter axis; and each 
row, as indicated by the number to the left of it, contains systems with a constant 
number of lines running in the direction of the longer axis. A nodal ellipse 
counts as two in the latter case. It will be seen that two nodal classes occur 
alternately in the first row and again in the third row, and the other two classes 
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Normal nodal systems of free elliptical plates. 
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in the second row and again in the fourth row. ‘Three spaces at the top left-hand 


corner must remain vacant. ‘The diagram can be continued indefinitely down- 
wards and to the right. 


§4. THE NATURAL FREQUENCIES 

Table 1, which is in two parts, has been arranged so as to correspond with the 
nodal systems of Figures 1 and 2 respectively. The frequencies are expressed 
in terms of the gravest tone taken as unity. Particulars regarding the dimensions 
and actual lowest frequencies in plates used are also included. ‘The numbers in 
italics are approximate, having been estimated by extrapolation from curves 
derived from the earlier frequencies of any given row or column, That the 
frequencies are in general agreement with Chladni’s approximate results can be 
seen by comparing Table 1 with Table 2, which has been compiled from his 
data concerning plates of the same shape as those used in the present 
experiments. 


Table 1. Relative Frequencies of Normal Vibrating 
Modes of Free Elliptical Plates 


2/1 Ellipse (see Figure 1) 


0 1 2 3 a 5 6 
0 is 2°58 4-7 Hes} 10 
1 ee Sel 5:68 8:29 11 
2 4-25 6°57 9°43 12-6 
3 10:6 14 
4 V7 fap 


* Actual frequency of brass plate, major axis 12°68 cm., minor axis 6:39 cm., 
thickness 1:62 mm., was 438 c/s. 


5/4 Ellipse (see Figure 2) 


0 1 2 3 ot 5 6 7 
0 aly Boss 4:28 6:66 9°39 13 
1 1:07 DY) 4-34 6-8 9-6 
2 2-03 399 6-71 10 
3 4-42 TAL 10-7 14 
4 9-01 12 
5 14 


* Actual frequency of brass plate, major axis 12-60 cm., minor axis 10:18 cm., 
thickness 1:62 mm., was 414 c/s. 


Table 2. Chladni’s Vibration Frequencies (approximate)* 
2/1 Ellipse 5/4 Ellipse 


0 1 Dee 4 ® 2764 1 2:55 455 09 
1 (6a 35 6 8s 151 2-5 4:5 7 9-2 
mmis:6° 66° 9-4 13 - 17 +20 nt oak O28 411 16 18 

mel 10. 14 918 20 DS sla ae 19 

a1. 19 24 89 12 

5 14 


* Converted from chromatic scale notation. 
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§5. FURTHER REMARKS 


(1) Comparison with Circular Plate | 

On a circular plate, which is the limiting shape of elliptical plates, nodal — 
ellipses become circular and nodal hyperbolic lines become diameters. | 

The correspondence between elliptical and circular nodal systems may be | 
traced on Figure 2 by proceeding as follows: 

Study two rows at a time, the first with the second and the third with the 
fourth. Couple each design in the upper of the two rows with the one preceding 
it in the lower of the two rows. ‘Then each pair of nodal systems will be found to. 
correspond to a single nodal system on the circular plate (Waller 1938). 

For example, 0)|2 in the first row and 1|1 in the second row both correspond 
to the two-nodal diameter system of the circular plate, 0|3 and 1|2 to three diameters, 
and so on. 

Similarly pairs of nodal systems can be picked out in the third and fourth 
rows, omitting the nodal ellipse which occupies the first place and which corre- 
sponds to the nodal circle. ‘Thus the 2]1 and 3|0 systems both correspond with 
a nodal circle with one diameter, and so on. 

The ratios of the frequencies of paired nodal systems can be studied in Table 1, 
noting how much nearer they are to unity for the wider than for the narrower 
plate, 1:07 for example as compared with 1-77 for the first pair. In the limiting 
case of the circle the ratio becomes equal to unity. 


(11) Comparison with Rectangular Plates and with Bar 

As previously noted, the nodal systems of elliptical and rectangular plates 
are arranged in the same way. ‘The nodal designs evidently become more nearly 
alike as the plates become narrower. On very narrow plates the only nodal 
systems met with in practice are those which can be placed in the first row of 
Figure 1. 

These approximate to the nodal systems of free bars where the lines run at 
right-angles to the length and are divisible into two classes according as the centre 
is antinodal or nodal. 


(11) Combined Modes of Vibration 

We may ask why among the photographic records of Figure 2 the 1|4 and 23 
designs are so distorted. ‘Turning to the second part, Table 1, it is found that 
the corresponding vibration frequencies are nearly equal, 6-8 and 6-71, respectively. 
The fact that the two modes have attempted to combine, somewhat imperfectly, 
indicates that the plate is not quite uniform and explains why the designs are so 
irregular. 

Incidentally this provides by far the most sensitive means of detecting a want 
of uniformity in a plate. But unfortunately it is unlikely to be of any practical 
use since combination occurs so rarely and cannot be predicted. 

We will conclude by mentioning one important difference between the nodal 
designs of rectangular and elliptical plates. 

On rectangular plates a trace of a second vibrating mode can sometimes be 
detected in a given nodal design, which always belongs to the same class of 
symmetry and which does not impair the symmetry of the design. Such 
combined vibrations have not been noticed on elliptical plates. This problem 
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of combined vibrations, due presumably to anticlastic coupling, is a difficult 
one. It can probably be studied most fruitfully by further experiments on plates. 
of different shapes. 
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Preparation of Alkali Metals in Glass 


When studying the electrical conductivity of the alkali metals it is necessary, owing to 
their great chemical activity, to prepare specimens for experiment in sealed glass capillary 
moulds. It is usual to make these of common ‘ soft’ glass so that platinum electrodes 
(which do not react with the metals or become oxidized during glass-blowing) may be sealed 
in. In recent work it has been necessary to use two special glasses for the preparation of 
potassium and lithium specimens respectively; we wish here to report our experiences 
with them. 

The electrical resistance of ‘ pure’ potassium specimens exhibited two distinct ‘ kinks ” 
in the low temperature region between 20° and 4° K. (MacDonald and Mendelssohn 1948, 
1950). To determine whether this was an intrinsic propeity of the ideally pure metal or 
due in some way to a small impurity content several distillations were carried out yielding 
purer samples. A limit, however, as judged by the residual resistance at very low temper- 
atures, was reached quite soon. It was presumed that this was due to impurity (sodium 
probably) being collected from the glass distillation container; in order to minimize this, 
‘Pyrex’ glass was used for making the container because of its lower sodium content. 
To effect further improvement a special glass (C 98) has been produced in which the only 
alkali ion present is potassium. This glass, of composition and properties shown in the 
accompanying Table, works well in a flame and can be drawn and blown easily. It is very 


Expansion coefficient 


=, B.0, Al,O, B30: -Ne,0°° K,0  —-Ei,O (0-400° c.) 
C98 64 2 = 18 = 16 nis 9-6 x 10-8 
C 51/65 69:2 = == 18-2 : = 12-6 10-3 x 10-8 
C 10 De 96-8 ~ 29:6 1350 4-4 1-9 oe 9-0 x 10-8 


Components are given in percentages by weight. 


resistant to devitrification (none appearing in five days at any temperature between 675° and 
1,000° c.) and may be sealed directly to ordinary soda-lime-silicate glasses of about the same 
expansion coefficient. When this glass was used for the distillation vessel a useful improve- 
ment in residual resistance was immediately observed and, furthermore, the ‘ kinks ’ were 
now found to have disappeared, suggesting very strongly that the resistance phenomenon 
was due fundamentally to small amounts of sodium impurity, although the fundamental 
physical source of the phenomenon is still obscure. 
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When attempting to work with lithium metal a serious obstacle is met with since the 
‘molten metal reacts very strongly with ordinary glass, forming almost instantaneously a 
black compound on the surface of the glass; the action generally proceeds so rapidly as to 
cause fracture of the glass itself. It has been suggested (see, for example, Justi 1948) that 
the high melting temperature of lithium is primarily responsible for this behaviour; since, 
however, sodium can readily be heated in glass to temperatures well above this value without 
appreciable reaction, this hypothesis seems untenable. 

It was thought at first that the process might be one of replacement of the sodium in the 


-glass by lithium, and, therefore, a special glass was made up which contained lithium as the 


only alkali ion. The well-known Lindemann glass is of this type (containing only lithium, 
beryllium and boric oxides), but it is difficult to work into the apparatus required by the 
experiments. The simple lithium-barium-silicate glass (C 51/65) of composition shown in 
the Table can be drawn into tubing (including capillary) by a skilled glass-worker, and 
although the glass has a low softening temperature and must be worked with care, the tube 
can be sealed directly to ordinary soda-lime-silicate glasses of similar expansion coefficient. 
“Thin platinum electrodes can be sealed into the glass without cracking. 

Unfortunately, tests soon showed that this glass was no more resistant to lithium than 
ordinary ‘ soft soda’ glass. Recalling, then, that fused silica is also attacked in the same way 
by lithium, it was thought that the silica in the glass was probably the responsible con- 
‘stituent. An experiment was therefore made using tubing of the type employed in the 
construction of sodium discharge lamps; the tubing has an inner lining of a borate glass 
which is almost free from silica, and is known to be resistant to the attack of both sodium and 
potassium vapours. The composition of the glass (C10) is also given in the Table. C10 


has proved sufficiently resistant to permit lithium to be cast successfully in it if the metal is 


‘solidified fairly rapidly (say 10 to 20 seconds) after running into the mould; thereafter a 
blackening does begin to occur. It is hoped in the future to make experiments with other 
glasses in which the silica content will be reduced to zero and the Na,O and K,O replaced 
iby, 11,0: 


Clarendon Laboratory, Oxford. D. K. C. MacDonaLp. 


British Thomson-Houston Co. Ltd., J. E. STANWORTH. 
Research Laboratory, Rugby. 
4th April 1950. 


Justi, E., 1948, Ann. Phys., Lpz., 3, 183. 
-MacDona_p, D. K. C., and MENDELSSOHN, K., 1948, Nature, Lond., 161, 972; 1950, Proc. Roy. 
Soc. A (in the press). 


The Dissipation of Energy by a Pendulum Oscillating in 
Air at Low Pressures 


The energy dissipation of a seconds pendulum has been measured by Atkinson (1938), 
who described experiments with different pendulums by means of which losses due to 
suspension, rod and bob were calculated. Such losses related to pendulums swinging in 
the open air, and these were found to be higher when the pendulums were inside cases. 
The loss of energy is due to four causes : (a) air resistance to the bob, (d) air resistance 
to the rod, (c) bending of the suspension spring, and (d) movement of the support. The 
combined loss due to causes (c) and (d) can be determined directly from measurements 
made with the pendulum in a vacuum. 

The air pressure in the cases containing the ‘ Shortt clocks’ is normally reduced to a 
value 30-20 mm. Hg, and it was decided to measure the losses of one of these pendulums 
at lower pressures and, by extrapolation to zero pressure, determine the losses due to the 
flexure of the suspension spring and the movement of the support. 
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The seconds pendulum of the standard Shortt clock seconds type, used for the tests, 
was carried by a spring ground from a solid strip of steel 2 mm. thick. The upper 
thick end of the strip was gripped between the two halves of a cylindrical trunnion while 
the lower thick end carried a steel pin which engaged with two hooks formed at the top end 
of the 8 mm. diameter invar pendulum rod. ‘The bob, made of type metal, was a cylinder 
9 cm. in diameter, 10 cm. long, weighing 6:5 kilograms. It was supported a little below 
its centre of gravity by means of an invar sleeve—pinned to the pendulum rod—carrying 
a brass collar of such a length that its upward expansion due to temperature compensated 
the corresponding downward expansion of the suspension spring and pendulum rod. 
‘The case was a copper tube of internal diameter 21-5 cm., fitted with heavy bronze flanges 
at the top and bottom carrying lugs by means of which it was solidly bolted to the wall. 
The trunnion, holding the top of the suspension, was supported by a four-legged bronze 
casting which stood on the inner edge of the top flange. This casting, which also carried 
the maintaining mechanism, was covered by a glass bell-jar, sealed to the outer edge of the 
flange. ‘The base of the case was closed by a thick glass disc sealed to the bottom flange. 
When the pendulum was swinging, the underside of the cylindrical bob was 8:2 cm. from 
this plate, but the rod was extended downwards to carry a silver ‘ beat plate ’ which moved 
to and fro just clear of the glass. This beat plate was engraved with a series of cross-lines 
spaced 1:67 mm. apart which was the appropriate interval to represent 5’ of arc. The 
movement of these lines, as the pendulum oscillated, was observed by means of a microscope. 
When illuminated by a suitably placed lamp, one of the plate lines came into view at the 
end of a swing, slowed to a standstill, and then swung ‘ out of focus’ again, followed by a 
similar movement of another line entering from the opposite side. It was quite a simple 
matter to read off these ‘ standstill positions’ on the eyepiece scale at any particular time 
or, alternatively, to observe the times when the lines ‘ stopped ’ at particular scale divisions, 
or at the same division. 

The rate of energy loss of a swinging pendulum is Jwdw/dt, where I is the moment of 
inertia of the pendulum and w the angular velocity when passing through its equilibrium 
position. ‘The value of J for the pendulum used was 6-852 x 10’ G.c.s. units and, repre- 
senting the total arc of swing by A, the rate of energy loss in ergs per sec. was 14:28A dA/dt. 
To determine the relation between the rate of energy loss and the gas pressure it was only 
necessary to measure dA/dt at 100’ total arc swing at various pressures. 
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The air pressure in the case having been adjusted to the required value, the arc of the 
pendulum was raised to 110’ or more by utilizing the normal impelling mechanism; the 
impulsing was then cut off. ‘The times were recorded at which the total arc became 105’ 
and 95’, respectively, and dA/dt determined. ‘The gas pressure was measured by means of 
an oil manometer at the higher pressures and by a McLeod gauge at the lower values. 
The results obtained are shown in the accompanying diagram. 

They indicate that the energy loss due to the flexure of the particular spring used, 
together with the losses due to the movement of the suspension brackets, i.e. of the head of 
the case, is less than 0-27 erg per sec., and probably reaches a limiting value of 0-20 erg persec., 
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that the damping at 3 x 10-4 mm. Hg is only about one-sixth the damping at 10-1 mm. Hg, 
and that gas viscosity damping still persists at a pressure of 10~4mm. Hg, which agrees 
with the fact that the quartz fibre vacuum gauge is operative at this pressure. The limiting 
value of 0-20 erg per sec., which represents the suspension losses, is considerably less than 
that previously estimated by the author—a value criticized by Atkinson as being too high— 
but the earlier figure was based on tests which were confined to pressures above 1 mm. Hg. 
The present tests also confirm the statement by Loomis, quoted by Atkinson, that the 
decrement of a half-second pendulum was reduced to one-half by decreasing the gas pressure 
from 0:025 mm. to 0:001 mm. He. 


Department of Physics, W. H. SHortr. 
University College, Exeter. 
5th April 1950. 


ATKINSON, E. C., 1938, Proc. Phys. Soc., 50, 721. 
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REVIEWS OF BOOKS 


Higher Physics, by E. NIGHTINGALE. Pp xvi+808. (London: G. Bell and 
Sons Ltd., 1948.) 27s. 6d. Also in separate parts: Pt. 1. Mechanics and 
Properties of Matter, 7s. 6d.; Pt. Il. Heat, 7s. 6d.; Pt. IL. Light and Sound, 
10s.; Pt. IV. Electricity, 12s. 6d. 


The clarity of style which has characterized Mr. Nightingale’s many school texts in the” 
past is 1n evidence in his latest effort. This book is of intermediate standard and is primarily 
intended for sixth form scholars. 

The wide field of subject-matter which is covered within the confines of this single 
volume has necessitated considerable ‘ condensing ’ in certain sections, but on the whole this 
process has been carried out quite efficiently. However, it is inevitable that some sub- 
divisions become unduly compressed when room is found for such topics as the theory of 
flight, the quantum theory of spectra, etc. 

An undesirable effect of the lag, so prevalent in the post-war era, between the preparation 
of an MS. and its publication is exemplified in the present instance by the author’s failure 
to record the recent official acceptance (October 1946) of the absolute system of electrical 
units as the practical system, to the exclusion of the International Scale. Also no mention 
is made, rather surprisingly, of the M.K.s. system of units. 

The author makes frequent use of the methods of dimensions in the various sections of 
the book and the use of graphs in exhibiting experimental results is also stressed in the text. 

The errors noted are remarkably few, but one which should be noted is the mention, 
on page 633, of the use of an iron former for the suspended coil of a galvanometer. 

Line diagrams which are capable of reproduction by the student are freely employed 
thoughout the text, and another feature which commends the book as one to be treasured by 
the scholarship candidate is the large selection of examination questions at the end of each 
chapter. R. W. B.S. 


. 


Terrestrial Radio Waves: Theory of Propagation, by H. BRemMeR. Pp. x+343. 
(New York, Amsterdam, London, Brussels: Elsevier Publishing Co., 
Inc.; British Agents: Cleaver-Hume Press, 1949). 


The author asks us in his Preface to bear in mind that his aim is to describe the 
mathematical-physical methods for the computation of transmitter fields. The statement 
of the problem is presented briefly in the opening pages and looks deceptively simple : it is 
required to compute the radiation field of a Hertzian dipole above a spherical earth. The 
influence of the ionosphere and the refracting lower atmosphere is deferred until Part II of 
the book. It takes some hundred odd pages of concisely written mathematical analysis to 
arrive at the final formulae for the first part of the problem in a form suitable for numerical 
work, and the reader will readily agree with Nicholson that the problem was one of the most 
difficult at one time facing the theorist. The author and his colleague Professor van der Pol 
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‘Were foremost in reducing its solution to a practical form for numerical work ; their original 

research papers, published in the Philosophical Magazine in 1937-1939, have long been 
recognized as the ‘‘ textbook’ of workers in this field. Indeed, much of Part I will be 
familiar to those who have read these papers. 

The apparent ease with which the author systematically takes us through pages of 
complicated mathematical analysis could only come after long experience and familiarity 
with his subject; his emphasis on the physical interpretation at each stage is an outstanding 
example of the presentation of mathematical physics. Even so, it is to the experienced 
theoretician that his descriptive passages are directed, and it is unlikely that the average radio 
engineer has the necessary mathematical equipment for their appreciation. Nevertheless 
the debt owed by the radio engineer, as shown by his extensive use of the results of this work 
in the past, was aptly expressed by T’. L. Eckersley in his concluding remarks to the Radio- 
location Convention in 1946 when he said that “‘ there is nothing so practical as a good 
theory’. 

Without much preliminary, a beginning is made to the solution of the point source and 
spherical earth problem. Starting with the classical solution, in the form of a series of 
zonal harmonics, the author takes us through the various transformations and approximations 
with systematic thoroughness, eventually leading to the familiar residue series. An instruc- 
tive account of the theories of steepest descent and stationary phase is of general interest in 
the relationship of ray to wave theory, and is more clearly discussed here than in most books 
on the subject. The author does well in giving prominence to the physical significance of 
geometric-optical approximations and of Watson’s transformation, which might easily have 
remained as so much mathematics. 

The concluding chapter of Part I will be particularly welcomed by the radio engineer 
whose interest is in the calculation of field strengths. The specific examples chosen and 
worked through in detail will be of assistance in using the formulae developed earlier in the 
book and now collected together in a convenient form for computation. Curves covering a 
wide range of frequencies are also included. 

Part II of the book is devoted to a discussion of the influence of the ionosphere and 
refraction by the troposphere. With the simplifying assumption of a stratified atmosphere 
whose refractive index varies only with height, a comprehensive discussion is given of both 
the residue series and geometric-optical solutions. The reader will find an enlightening 
account in Chapter VII of the method by which the spherical earth problem can be replaced 
by one with a “‘ flat’’ earth and a modified index of refraction. Closely allied to this is 
Eckersley’s method of taking refraction into account to a first approximation by using a 
“ fictitious ”’ radius of earth which is usually larger than the actual radius. During the war 
years this theory has been supplemented by the more refined mode theory of super-refraction 
which explains the abnormally long ranges experienced by microwave radar during favour- 
able meteorological conditions. "The complexities of the eigenvalue problem involved in 
the evaluation of height functions receives detailed attention ; the discussion of the W.K.B. 
and the more accurate Krammer phase integral approximation is of general interest to the 
mathematical physicist. 

This brief outline of the book cannot possibly do justification to its wealth of detail. It 
is without doubt a book for the propagation specialist who requires a fundamental under- 
standing of his subject. There is, however, much of general interest for the theoretician. 

/ W. WALKINSHAW. 


Oscillations of the Earth’s Atmosphere, by M. V. WiLKes. Pp. 1x+74. Cam- 
bridge Monographs on Physics. (Cambridge: University Press, 19449). 
12s. 6d. net. 


This excellent little volume by the Director of the University Mathematical Laboratory, 
Cambridge, deals with the oscillations produced daily in the atmosphere by the sun and the 
moon, as revealed mainly by the solar and lunar daily variations of barometric pressure. ‘The 
subject may be said to have started with Torricelli’s invention of the barometer and Newton’s 
gravitational theory of the tides, greatly developed by Laplace, who in his last years devoted 
much attention to atmospheric oscillations. Kelvin showed the importance of resonance 
in greatly magnifying the semidiurnal atmospheric oscillations due to the sun’s thermal and 
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tidal action, so that the lunar tide in the barometer, unlike that in the sea, is relatively 
inconspicuous. The present data on the subject are described, and also Laplace’s tidal theory. 
The last and most novel part of the book is Chapter V which gives the modern theory of 
these oscillations and their resonance, a theory begun by Taylor and Peteris and since then 
much developed by the author and Dr. K. Weekes. The book should interest many 
physicists, mathematicians and meteorologists and stimulate further developments of which 
the subject still stands in much need. S.C. 


Sound (A University Text-book of Physics, Volume II), by J. H. PoyNTING and 
J. J. THomson, revised by W. S. Tucker. Pp. viii+251. ‘Tenth Edition. 
(London: Charles Griffin, 1949.) 20s. 


During the last 50 years nine editions of this well-known textbook have been published. 
It is not surprising that in this present edition the reviser has found it necessary to write 
effectively anew book. The text has been rearranged and the theory simplified by the use 
of calculus methods. The reviser claims that the book has been written to remedy a 
deficiency in present University courses, where ‘ Sound’ “‘ is frequently relegated to the 
end, and is largely crowded out’’, and has been planned for students taking honours courses 
in branches other than physics, viz. those who have chosen physics as a subsidiary subject 
for their degrees. 

The text mostly divides into two parts—pure and applied acoustics. The former 
follows normal procedure and requires little comment except that an improvement could be 
made in the definition of terms introduced and in the physical—mathematical analysis 
relation. In the applied acoustics part the question of relevancy of subject matter is much a 
question of opinion but it would appear that the text would gain in clarity and cohesion by 
considerable pruning and rearrangement. A considerable proportion of the book, about 
34 pages, is devoted to the hot wire microphone and its application to sound location of 
aircraft and guns, and a further 10 pages to the acoustical design of aircraft sound locators. 
A discussion of the siren takes thirteen pages. Although aircraft noise and its measurement 
is considered in eight and a half pages there is little if anything on other noises and their 
measurement and suppression. Certain subjects are introduced in an unusual way, e.g. 
the Doppler effect through aircraft sounds without further mention of application, and the 
Whispering Gallery effect through observations on large sound reflectors for the detection of 
approaching aircraft. Some of the treatment appears to be unnecessarily dated, e.g. the 
recorded frequency spectrum of gramophone recordings is quoted as 50—5,000 c/s. and the 
last reference given is 1926. Nomention is made of the recent improvements in wire and tape 
magnetic recording. Readers may argue however that such details are unnecessary in a 
textbook of this standard since they do not add to the fundamental concepts. 

An erroneous impression may be given by certain statements such as :; on p. 93, the uses 
of the tuning fork as a present standard of frequency for radio and clock control; on p. 242, 
the highly absorbent walls of broadcast studios; and on p. 134, the description of the 
mechanics of the electronic musical instruments. There seems to be some misunderstanding 
on the relative positions of the sound and vision gates on film projectors. 

On musical questions, scale, consonance, etc., tend to be considered as a result of numerical 
arrangement. Numerical approach has interest but is mostly sterile. 

The last chapter in the book, viz. ‘‘ The Acoustics of the Concert Hall ”’ could be better 
titled. Only a quarter of the printed matter refers to this question. 

In a future edition an improved text could be obtained by a rearrangement of the 
material, with the deletion or considerable reduction of certain sections and the incorporation 
of further diagrams and photographs on the pure acoustical side. The addition of some 
thought-provoking problems at the end of chapters is also suggested. H. D. P. 
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High Vacuum Technology, by A. S. D. Barrerr. 


Abridged version of the 23rd Anniversary Address given to the King’s College 
Engineering Society at the Institution of Mechanical Engineers on 10th March 1949. 
(London: W. Edwards & Co. Ltd.) 
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ABSTRACTS FOR SECTION A 


The Schumann—Runge O, Emission Bands in the Region 3100 a.-2500a., by 
Nie Wo Pens. 4 


ABSTRACT. Rotational and vibrational analyses of the emission Schumann—Runge O, 
bands in the region 3100 A. to 2500 A. are presented. ‘Tables of wave-numbers are given for 
the (1-12), 11), (4511), (0; 10), 10); (5.9); 2,9), (Gy 8)e 28) (2597) eibanicissaeeeeens 
origins of the bands, the rotational differences, and constants for the levels v’=7 to 11 are 
given. ‘Three bands, found in a high voltage arc in oxygen between platinum electrodes, 
are attributed to PtO. 


New O,+ Second Negative Bands: A Note on O, and Ou Emission Spectra, by 
M. W. Feast. 


ABSTRACT. The spectrum of an electrodeless discharge in pure oxygen has been 
studied. New bands of the O,+ Second Negative system have been found and the intensity 
distribution in this system is discussed; subsidiary parabolae as well as the main Franck— 
Condon parabola are obtained as is predicted by the theory. Emission bands reported by | 
Johnson as due to Og; are due to O,+ and O11. Some abnormal intensities of Ort lines | 
excited in the electrodeless discharge are noted. The electrostatic and electromagnetic 
types of excitation are briefly compared. 


| 


The Spectra Emitted by the High Voltage Arc in Nitrogen, Hydrogen, Nitrogen~ 
Hydrogen Mixtures and Ammonia, by M. W. Feast. 


ABSTRACT. The spectra emitted by the high voltage arc in nitrogen, hydrogen, 
hydrogen-nitrogen mixtures and ammonia-have been studied. 

The appearance of the N, Second Positive bands at various nitrogen gas pressures is 
discussed and the divergence of rotational and vibrational temperatures is noted. Traces 
of oxygen greatly reduce the strength of the N.* bands emitted by the arc in nitrogen. 
In N,-H; mixtures the NH 3360 a. band is found to be emitted strongly only when hydrogen 


is present at a very low partial pressure. An attempt is made to explain the various 
observations. 


Rotational Analysis of the (1,0) Band N, First Positive System, by M. W. Frast. 


ABSTRACT. The (1, 0) band of the N, First Positive system has been photographed at 
a dispersion of about 1-2 cm-!/mm. The wave numbers and quantum numbers of the 
lines, and the rotational term differences in the upper and lower states are given. The 
spin tripling in the *% state and the A doubling in the °I] state are evaluated. 
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An Absorption Comparison. of the B-Particle Spectra of *%' AcC" (allowed), !8RaE 
(second forbidden) and 3°5 yr.-*3{ TI (third forbidden), by H. D. Evans. 


ABSTRACT. Experiments have been described in which a direct comparison is made by 
an absorption method of the f-particle spectra of actinium C”, radium E and 3:5 yr. thallium 
representing allowed, second forbidden and third forbidden transitions respectively. 

New values of 623-+-4 mg/cm? and 300+3 mg/cm? are proposed for the absorption 
limits in aluminium of the B particles emitted by actinium C” and thallium and it is shown. 
that, of the three elements, radium E emits most and actinium C” least low energy electrons. 

An attempt has been made to derive energy distributions from the absorption curves, but 
the results for radium E depart markedly from those obtained with magnetic spectrographs. 
The general findings, however, agree with those of other workers who have attempted the 
same procedure. 

The formula used in the derivation of the energy distribution includes an approximation 
which necessarily leads to an apparent deficiency of low energy particles, but the accuracy of 
the data is not sufficient to permit the use of a better approximation. Strictly for purposes of 
comparison, however, the method has obvious applications. 
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Excited Electronic Levels in Conjugated Molecules—V : A Valence Bond Estimation 
of Energy Levels in Aromatic Hydrocarbon Molecules, by G. R. BaLDock. 
ABSTRACT. This note is concerned with the calculation of the energy of the ground 
state and the lowest excited state for hydrocarbons of the condensed ring type. The differ- 
ence between these two energies is the energy of the first electronic transition. "The normal 
valence bond method is used, taking only the Kekulé structures into consideration. The 
results are compared with the experimental absorption wavelengths, and are in good 
agreement, giving a value of the exchange integral « between —1-8 and —2 electron volts. 


Relation between Bond Structure and the Longest Ultra-violet Absorption Band of 


Polycyclic Aromatic Hydrocarbons, by J. W. Coox, R. ScHOENTAL and 
Pi. ¥. SCOTT: 


ABSTRACT. A comparison has been made between the values expressed in wave 
numbers of the longest ultra-violet absorption bands of a series of polycyclic aromatic 
hydrocarbons and the minimum number of quinonoid rings contained in their molecules. 


It is found that the hydrocarbons fall into series of groups related to this number of quinonoid 
rings. 


The Scattering of High Energy Charged Particles by Thin Foils of Matter, 
by S. T. BuTier. 


ABSTRACT. The distribution of high energy charged particles scattered by a thin foil 
has been investigated and compared with the formula for single scattering. It is shown that 
at sufficiently large angles the distribution approaches that for single scattering and can be 
expressed in the form of a rapidly convergent series of which the first term is identical with 
the single scattering formula. The remaining terms of the series form an accurate expression 
for the deviation from true single scattering. The simplified case of projected scattering has 
also been considered. 


Penetrating Showers at High Altitude, by M. G. Noou anp S. R. Happara. 


ABSTRACT. A study has been made of cloud chamber photographs of penetrating 
cosmic-ray showers, originating in lead placed over the cloud chamber and in a lead plate 
within the cloud chamber, obtained at altitude 3,572 m. in a magnetic field of about 1,000 
gauss. 

The use of well-separated successive lead layers allowed the association of penetrating 
events to be demonstrated, and there is some evidence of a general reduction of energy 
involved in events in the lower lead layers as compared with those in the top layer. Mesons, 
protons and particles of higher charge were identified in the penetrating showers, and a 
mean range for non-coulomb scattering of shower particles is given. Evidence is also 
given of the apparent starting and stopping of charged particles in the lead plate at the centre 
of the cloud chamber. 


On the Stresses and Energies associated with Inter-Crystalline Boundaries, by 
J. H. van per MErwE. 


ABSTRACT. Models, largely base don the assumptions introduced by Peierls and 
Nabarro in dealing with a single dislocation, are used in calculations on three types of 
intercrystalline boundaries, namely, (I) a boundary due to a difference of atomic spacing, 
(II) a twist boundary, and (III) a symmetrical tilt boundary. With these models the 
resolution of the boundary into a sequence of dislocations is a natural consequence of the 
analysis, which also yields the expressions for the stresses, atomic displacements and energies 
associated with the boundary, as functions of the angle of tilt, etc. By allowing the distance 
between dislocations to tend to infinity, these expressions reduce to the corresponding 
ones for single dislocations. The outstanding feature of the interfacial energy is that it 
increases initially very rapidly with the angle of tilt, etc. An application of the results to 
the theory of oriented overgrowths, developed by Frank and the present author, is described. 
The validity of the assumptions and approximations involved and the advantages of the 
treatment as compared with those of other writers are discussed. 
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The Theory of the Propagation of First and Second Sound in Helium II.—Energy 
Theorems and Irreversible Processes, by R. B. DINGLE, 


ABSTRACT. The first part of this paper contains a discussion of the conditions deter-_ | 


mining whether sound waves in helium II are propagated isothermally or adiabatically. In 
the remainder of this part an expression is found for the energy flux from a closed region, and 
the impedance concept is applied to waves of second sound. 

The second part of the paper is devoted to the study of the effects of irreversible processes 
on the propagation of waves of first and second sound. ‘The irreversible processes discussed 
are those due to viscosity and thermal conduction, both in an extended fluid and in a narrow 
tube. 


The Angular Distribution of Synchrotron Target Radiation: A Preliminary 
Experimental Study, by J. D. Lawson. 


ABSTRACT. The angular distribution of radiation from a synchrotron has been measured 
with an ionization chamber for a number of different target materials and thicknesses, and 
the results have been compared with a simple theory in which the radiation from a parallel 
beam of electrons striking a thin plate has been calculated. Detailed comparison with theory 
is not possible, because the precise way in which the electrons hit the target cannot be 
determined, nevertheless qualitative explanations have been found for most of the effects 
observed. The most important of these is the relative insensitivity of forward intensity and 
beam width to target size and material. : / 
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